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A calibrated spectrometer can indirectly determine a physical or chemical property of a sample based upon spectral re- 
sponses measured by the spectrometer with respect to the particular sample. This invention is directed to a method for calibrating 
or recalibrating a second spectrometer in light of a first spectrometer, or itself, respectively. The calibration employs a unique se- 
lection and manipulation of spectral data obtained from both the first and the second instrument. The recalibration employs a un- 
ique selection and manipulation of spectral data from the same first instrument, that is obtained both before and after the need 
for recalibration arises. Instead of modifiying the respective responses of the first and second instrument, or the first instrument 
before and after the need for calibration arises, this invention modifies the calibration equation of the second, or recalibrated in- 
strument, to yield consistent results to those obtained by the first instrument, or the first instrument before it goes out of calibra- 
tion. A calibration equation is an equation which transforms spectral data of a particular sample at a variety of wavelengths to a 
calculated value for a chemical or physical property. Generally, the form of such calibration equations is that of a linear combina- 
tion of absorbances or mathematical transforms of absorbances measured by the first and second instrument for each sample. 
Surprisingly, and confirmed statistically, the accuracy of the calibrated or recalibrated instrument is maintained and in some in- 
stances improved. 
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2 . V 

* SPECTROSCOPIC INSTRUMENT CALIBRATION 

5 BACKGROUND OF THE INVENTION 

6 I. Field of the Invention 

7 This invention relates to spectrophctcmeters capable of measuring 

8 electromagnetic radiation, such as infrared, near-infrared, ultraviolet, visible, Ramon 
9 

1Q spectra, Rayleigh scattering, and the like. More specifically, this invention relates 

^ to calibration techniques useful for predicting chemical or physical properties 

12 based on observed spectral data, preferably data obtained by means of a 

13 spectrophotometer. Still more preferably, this invention embodies near-infrared 

14 calibration equations suitable for making predictions or estimates of physical or 

15 

16 chemical properties such as octane. 

17 li. Prior Art Summary 

18 In J. Sci. Food Agric, 34(1983), p. 1441-1443, B. Q. Osborne directly 

19 transfers a calibration equation from one instrument to another. 

20 'Use of Indirect Multivariate Calibration for Qualify Control of Agricultural 

21 Products by Near-Infrared Spectroscopy,' G. Puchwein and A. Ebefiuber 

22 

(Mikrochim. Acta (Wien), 1986 II, 43-51) teach a computer program for ceJttxating 

24 an instrument. The program picks relevant frequencies for spectral data for part 

25 of a set of samples with which to calibrate an instrument and then validates with 

26 respect to the remaining samples. 

27 'Near-Infrared Limitations to Silicon Photddetector Self-Calibration,' 

28 J. M. Palmer (SPIE Vol. 499, Optical Radiation Measurements (1984], p. 7-14) 

29 

30 teaches self-calibrating silicon photodiode detectors. 

31 "A Simple Internal Modulation Technique for the Spectral Calibration of 

32 Circular Variable Filter Spectrometers In the Near-Infrared,' Chandrasekhar, Ashok, 

33 Bhatt, and Manian (Infrared Physics, Volume 24, No. 6, 1984, p. 571-572) teach 

34 the wavelength calibration of telescopes attached to near-infrared instruments for 

35 astronomical studies using the frequency dependent signal of a low-pressure 

36 

37 mercury arc lamp. 

38 Mark and Workman (Spectroscopy, 311), 28) developed a computer 

39 algorithm to find isonumeric wavelengths. Isonumeric wavelengths are calibration 
40 
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wavelengths which yield the same absorbance readout on a master NIR and a 

5 slave NIR. These wavelengths therefore do not change substantially from one 

6 instrument to another. 

7 Heckman, Dtffee, and Mllhous (Anatytica Chemica Acta, 192(1 987), 197) 

8 teach a computer program that transfers calibration equations by computing slope 
9 

1Q and bias corrections for each wavelength in the equation. This method of slope 

n and bias corrections is, in fact, used in Example 3 of this application, in 

12 US 4,761 ,522 R. D. Rosenthaul claims a poly-tetrafluoroethylene standard for 

13 near-infrared and processes for making that standard. That standard, however, 

14 is not used in calibration transfer from one instrument to another. 

15 Multivariate Instrument Standardization, by Wang et al., Anal. Chem. 1991, 

16 

1? 63, 2750-2756, discloses four standardization methods for either correcting the 

18 primary calibration model for use on a secondary instrument, or correcting the 

1 9 spectra acquired on the second instrument to account for the response differences 

20 between instruments. It was discovered that the calibration model built from one 

21 batch might not be applicable to another batch. However, using mathematics of 

22 

multivariate calibration, four different approaches were followed to select the best 

23 

24 subset of calibration transfer samples available. Among the methods also 

25 discussed was methods for correcting the spectra by a process of linear transfer 

26 of one spectral response in its intensity from instrument to that of another. Not 

27 disclosed in this reference is the use of a calculated dependent variable 

28 corresponding to a physical property of interest being used as the key to 

29 

3q establishing coefficient for constants in the second calibration equation. 

31 B. G. Osborne, in an article entitled, Calibration of Instrument for 

32 Near-Infrared Spectroscopy, (48 Spectroscopy, Vol. 4 No. 4, p. 48-55). Osborne 

33 discloses in detail the steps for establishing simple correlations by multiple linear 

34 regression and also the use of derivatives of linear terms. Validation of calibration 

35 curves are discussed; specifically how particular wavelengths are chosen to 

36 

3? establish the terms in the multiple linear regression equations, along with the 

38 statistical significance of a particular calibration sample. This method can be 

39 utilized in this invention for the purpose of defining the calibration equations of the 

40 
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calibrated instrument for a range of values within which useful predictions can be 

5 made, i.e., the range of interest. 

6 Still other methods for determining calibration equations are those discussed 

7 in Analytical Chemistry Journal voi. 60, pages 1133-1202 (1988) by David M. 
Haaland and Edward V. Thomas. The principles of inverse least squares, partial 
least squares, classical least squares, and principal component regression analysis 

y 1 are disclosed. 

12 Uquid Absorption Standards For Ultraviolet, Visible, and Near-Infrared 

13 Spectrophotometry, R. G. Martinek, J. Amer. Med. Technol. July-August, 1978, on 
p. 210-216. Liquid absorbance standards are shown to have the following 
desirable characteristics: broad absorption peaks; wide wavelength range; high 

1 ? molar absorptivity; stability in solution; readily definable specifications of purity; and 

18 minimal spectral temperature coefficients. The purpose of using such standards 

19 is to permit calculation of absolute absorbance based upon the observed 

20 absorbance and the application of an instrument correction factor. 
The Use of Statistics In Calibrating and Validating Near Infrared On-Une 

Analyzers, Bruce Thompson, ISA, 1989 Paper #88-0116, uses near-infrared for 

24 both process control and statistical process control. 

25 'Vibration of Instruments for Near-Infrared Spectroscopy", B. G. Osborne 

26 (Spectroscopy, 4(4), p. 48-50) teaches the need for multiplicative scattering 
correction for solid samples due to particle-size differences in solids and computer 
transfer of complex equations. Osborne discusses procedures for developing 
calibration equations of particular relevance and usefulness in the present 

3 1 invention. The selection of appropriate wavelengths and samples are discussed 

32 and relevant articles cited further explaining the techniques employed. The transfer 

33 of calibration equations from one instrument to another are disclosed to require, 
in the case of most modern instruments, only a change in Intercept value, and 
computerized procedures allow for differences between spectral data generated 
by different instruments. These computerized procedures are not taught in detail, 

38 but are likely to involve transforming the spectral data of one instrument to be 

39 consistent with that of another instrument. This is to be distinguished from the 

40 
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method involved In the present Invention, which instead of transforming the 

4 

5 spectral data, the relevant coefficients in the calibration equation for the second 

6 instrument are changed from those found in the calibration equation of the first 

7 instrument in accordance with a unique procedure discussed in this specification. 

8 In the inventor's doctoral dissertation, A Chemometric Analysis of a 
' J Magnetic Water Treatment Device available from University Microfilms 

1 1 International, Dissertation Information Service, order number 8919928, there is a 

1 2 discussion of linear regression analysis (simple and multiple) for a linear regression 

1 3 model, and data transformations which can transform non-linear types of data into 

14 data which is appropriate for linear regression modeling. Definitions are shown 

15 which establish the degree of reliability that the value of one variabi (the 

16 

17 dependent variable) can be estimated on the basis of values for other variables 

18 (the independent variables). In general, the methods for determining regression 

19 constants (or coefficients) for the calibration equations discussed in the inventor's 

20 doctoral thesis are useful in establishing the calibration equation for the first or 

21 reference Instrument in the present invention. 

22 

23 A New Approach to Generating Transferable Calibrations for Quantitative 

24 Near Infrared Spectroscopy, H. Mark and J. Workman, Spectroscopy Volume 3, 

25 No. 11, p. 28-36, shows that instrument variability makes transfer of calibrations 

26 between instruments difficult However, by selecting the wavelengths and 

27 corresponding spectral properties for such wavelengths to meet certain criterion, 

28 many of the variations between instruments will not change significantly the 
29 

30 calibration equation from one instrument for that calibration equation which would 

31 work satisfactorily in another instrument Not disclosed is how to revise one 

32 calibration equation from one instrument to be suitable for use in another 

33 instrument where there is a significant change in the calibration equation suitable 

34 for the second instrument. 

35 Tobacco constituents to Tilting Filter Instruments, R. A. Heckman, 

36 

J. T. Diffee, L A. Milhous, Analytica Chemica Acta, 192(1987) p. 197-203, teach 

38 transferring a near-infrared monochromatic calibration for tobacco constituents to 

39 a titting-filter instrument in order to utilize the experimental calibration equation for 

40 
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4 a sensitive but not rugged laboratory instrument to analyze data in a rugged 

s production oriented instrument to predict properties such as moisture, glycerin, 

6 propylene glycol, nicotine, and reducing sugar. A program called MTRAN cannot 

7 be used without difficulty uniess care was exercised to Insure that the ti!ting=f3ter 
instrument is property configured with filters so as to accommodate the 

9 

1Q wavelengths in the primary or master calibration equation for the laboratory 

1 1 instrument Consequently, before the primary calibration can be transformed, the 

1 2 monochromatic data files are compressed to simulate a tHting-filter instrument The 

13 resulting data file is used to re-establish the calibration equation. ' With the 

14 foregoing changes, CALTRAN, a program formerly used only for calibration 

15 

transfer between filter instruments, could be used to effect transfer from a 

16 

17 near-infrared monochromatic instrument to a tilting filter instrument. 3 

18 Typically, the prior art seeks to make changes in the observed spectra (or * - 

19 mathematical transform thereof) produced in one instrument conform to those ^ 

20 produced by a calibrated reference instrument. Only after the spectra from an 

21 uncalibrated instrument is brought into coincidence with that spectra from a 

22 calibrated instrument did one seek to transfer a calibration equation from one: 

24 instrument to another. Our experience has been that this tends to - introduce^ ^ 

25 errors, uncertainty, and numerous numerical difficulties which reduce the~ease and 

26 precision that is made possible by the invention of this specification. 

27 U.S. 4,963,745 to Maggard entitled Ocfane Measuring Process and Device 

28 discloses how to initially establish a calibration equation for a reference instrument 
29 

Specifically disclosed is that spectral features, such as second derivatives, of 

31 infrared spectral data at methyne band frequencies, in the range of about 1 200 to 

32 1236 nanometers (nm), were highly correlatable to octane values for gasoline-fike 

33 substances. Numerous patents and articles are cited within the text of the 

34 . specification, as references cited, and as other publications, all of which including 
35 

those expressly cited are expressly incorporated herein by reference. 

36 

3? An example of a linear caBbration equation is one of the form 

33 Y e ' - A' + £, B, 1 Xje 1 where a Y c ' is the predicted value of some physical or chemical 

39 property, such as octane values for a gasoline, C, of a set called sample-c. As 
40 
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discussed in Example 7, 226 gasoline samples were measured to determine A, 

4 

5 and each B^' for j = 1 to 3, to obtain a multiple regression equation suitable for 

6 predicting pump octane values of gasoline samples, using a first instrument; 

7 instrument-!, the reference instrument. Pump octane is the average of ASTM 

8 methods D2699 and D2700, and is commonly known as (R+M)/2 octane. The 

g 

io inferred reproducibility of this test method is a 95% confidence limit of ±0.71 octane 

n numbers over the octane range of 84-94 pump octane. From the second 

12 instrument, instrument-k, while in transflectance mode, gasoline absorbances at 

13 1 196, 1220, and 1236 nanometers (nm) were used to determine the values for the 

1 4 second derivative of an absorbance spectra, the spectral features, X^", in the linear 

15 

calibration equation for inferring a pump octane. The standard error of prediction 

16 

17 ("SEP*) is the statistical parameter that equals the square root of the mean 

1B residual variance of unknown samples. 

19 

20 SUMMARY OF THE INVENTION 

21 A calibrated spectrometer can indirectly determine a physical or chemical 
22 

property of a sample based upon spectral responses measured by the 

24 spectrometer with respect to the particular sample. This invention is directed to 

25 a method for calibrating or recalibrating a first spectrometer in light of a second 

26 spectrometer, or itself, respectively. The calibration employs a unique selection 

27 and manipulation of spectral data obtained from both the first and the second 

28 instrument The recalibration of the first instrument employs a unique selection and 
29 

30 manipulation of spectral data from the same first instrument, that is obtained both 

31 before and after the need for recalibration of the first instrument arises. Instead of 

32 modifying the respective responses of the first and second instrument, or the first 

33 instrument before and after the need for calibration of the first instrument arises, 

34 this invention modifies the calibration equation of the second, or recalibrated 

35 

instrument, to yield consistent results to those obtained by the first instrument 

33 

37 before it goes out of calibration. A calibration equation is an equation which 

3g transforms spectra) data of a particular sample at a variety of wavelengths to a 

39 calculated value for a chemical or physical property. Generally , the form of such 

40 
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calibration equations is that of a finear combination of absorbances or 

5 mathematical transforms of absorbances measured by the first and second 

6 instrument for each sample. Surprisingly, and confirmed statistically, the accuracy 

7 of the calibrated or recalibrated instrument is maintained and in some Instances 
improved. 

Instruments available on the market today take raw spectral data and 

1 1 digitize or form an analog signal representative of it for use in making calculations. 

12 As an essential step in such calculations, a calibration equation in conjunction with 

13 said signal defines a dependent variable. This dependent variable is either equal 

14 to or indicative of a value for some physical or chemical property of the material 

15 

^ from which the raw spectral data was measured. The calibration equation 

17 appropriate to a particular instrument will vary depending upon changes in the 

18 spectrometer that produces the raw spectral data. To calibrate a spectrometer, 

19 one defines an appropriate calibration equation that in conjunction with a signal 

20 representative of raw spectra measured by the spectrometer for some material will 

21 determine for the value of Its dependent variable a value either equal to or 

22 

indicative of some physical or chemical property of the measured material. In 

24 practice, changes in such an equation relate directly to changes in constants found 

25 in such an equation. This invention is directed in part to how to determine these * 

26 constants for an uncalibrated instrument by recourse to a calibrated instrument 

27 and its calibration equation. Once these constants are determined, then a new or 

28 transferred calibration equation has been found that can be used in conjunction 
29 

3Q with a digitized signal or an analog signal representative of raw spectra measured 

31 by the previously uncalibrated instrument 

32 For example, if one were to determine a calibration equation, calibration 

33 equation-2, for a second instrument, instrument-2, by recourse to a calibrated first 

34 instrument, instrument- 1 , and its corresponding calibration equations, equation* 1 , 

35 

the steps in any available order for achieving this in accordance with this invention 

36 

37 comprise: 

33 (a) determining from calibration equation- 1 a value indicative of each 

39 dependent variable- 1 corresponding to each member of a calibration 

40 
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4 the steps in any available order for achieving this in accordance with this invention 

5 comprise: 

6 (a) determining from calibration equation- 1 a value indicative of each 

7 dependent variable- 1 corresponding to each member of a calibration 

8 set based on spectral data-1 for each member of said calibration set 
9 

measured by means of instrument- 1; and 

10 

1 1 (b) determining constants-2 in said calibration equation-2 so that the 

12 sum of absolute differences between: 

13 (1) said value indicative of each said dependent variable-2 for 

14 each member of said calibration set based upon calibration 
15 

calibration-2 and data-2 corresponding to each member; and 

16 

17 (2) that corresponding value of dependent variable-1 calculated 

18 for each said member based upon calibration equation-1 and 

19 data-1 corresponding to each member; 

20 is substantially minimized. 
21 

22 

23 BRIEF DESCRIPTION OF THE DRAWINGS 

24 Figure 1 is a schematic overview of this invention. 

25 Figure 2 is a flow diagram for Example 5. 

26 Figure 3 is a flow diagram for Example 6. 

27 

28 DETAILED DISCUSSION OF THE INVENTION 

29 

30 This invention solves the problem of how to utilize the results from 

31 calibrating one instrument for calibrating a different instrument. Alternatively, this 

32 invention can deal with recalibration of an instrument to insure consistent results 

33 after there has been some change in an instrument's spectral response, such as 
occur due to changes in light sources, for example due to replacing burned out 
bulbs, or due to new connections used to collect data from different sources, such 

36 

37 as. changes in fiber optics used to transmit light or electromagnetic radiation from 

38 one location to another, due to changes in path length, or due to changes in 

39 spectroscopic mode, such as from transmittance to transflectance. 

40 
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The invention broadly deals with calibrating a new instrument by recourse 

5 to a reference instrument, but also contemplated is a method for continuously 

6 monitoring an instrument to be sure that it has not changed so much so that a 
new calibration equation is required to obtain reliable predictions for some physical 
or chemical property of a material of interest. 

The invention contemplates methods for indirectly determining a variety of 
1 , properties-p such as cetane, octane, aromatic-naphthene-olefin content, 

12 viscosity-index, molecular weight, and many others, by means of spectroscopic 

13 data. Spectral data broadly is intended to include, but not be limited to, all 

1 4 intensity versus wavelength spectra that arise when light, broadly intended to cover 

15 

a wide range of frequencies such as radio frequency to ultra-violet interacts with 

16 

1? a material. Generally, any spectral features (as defined in this specification) 

18 derivable from spectral data, which can be correlated to predict values for a 

19 chemical or physical property determined directly or indirectly, are within the 

20 intended scope of this invention. 
Some benefits made available by this invention include: 

1 . Only a small number of samples (compared to the number used in 

24 the original calibration of the reference instrument) is required to 

25 determine the calibration equation for the uncalibrated instrument 

26 because the calibration equation of the reference instrument can be 

27 used. Example 7 shows that a small calibration sample consisting 
of only four members can be used. 

2. The calibration set, sample-c, does not have to contain any unstable 

31 exemplars that may have been included in the original calibration set 

32 of samples and preferably, in the case of mixtures, such as gasoline 

33 or other petroleum fuels, can contain some of the pure molecules or 
easily prepared mixtures thereof. 

The limitation on appropriate members of sample-c in determining a 
37 calibration equation appropriate to an uncalibrated instrument based upon its 
33 spectral responses depends in part upon the statistical treatment chosen for the 
39 data. For multiple linear regressions, the number of members of sample-c must 
40 
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be at least one more than the number of wavelengths (independent variables) used 
in the model, unless a slope and bias correction or other error reducing step is 

6 also used. For example, at least four samples must be used to transfer a 

7 calibration equation having three independent variables. 

8 We have found that in addition to a slope and bias adjustment, other error 
reducing steps such as least squares error reduction, generalized least sum, and 
perhaps other techniques can further reduce any remaining error in the resulting 

12 calibration equation-2. 

13 Preferably, for the reference set (sample-c) to be useful for transferring the 

1 4 statistical variance, co-variance, or function derived from variance from a calibrated 

1 5 instrument to an uncalibrated instrument, or to the same instrument at a later date, 

16 

17 the sample-c should disclose spectral features which change as the instrument 

18 changes. For purposes of this invention, sample-c should be selected so its 

19 absorbance shows changes at the wavelengths of interest when changes in 

20 instrument performance occur, but should remain constant (within instrumental 

21 limits) when the instrument is operated at standard or typical conditions. 
A specific application of one solution derived from this invention is the 

calibration of a near-infrared spectrophotometer by use of the results derived from 
2s a fully or partially calibrated reference near-infrared spectrophotometer. 

26 Throughout the specification, reference will be made to one or more 

27 spectrophotometers. The fully or partially calibrated spectrophotometer or other 

28 instrument will be referred to as the 'reference instrument" or "reference 
spectrophotometer", has its calibration equation fully or partially defined. The 
instrument to be calibrated with respect to the reference will be referred to as the 

32 "uncalibrated instrument" or "uncalibrated spectrophotometer". A reference 

33 instrument may become uncalibrated or go out of calibration as a result of a 

34 variety of changes either within the instrument itself or from connections between 

35 the instrument and the source from which it derives its measurements or spectral 

36 

data. 

37 

33 This invention enables an uncalibrated spectrophotometer to be calibrated 

39 by reference to a reference spectrophotometer by means of a far smaller 
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4 calibration set of samples (the calibration set of samples) than was required to 

5 fully define a calibration equation for the reference spectrophotometer to begin 

6 with. The calibration set of samples is here referred to as set of sample-c, or 

7 simply, sample-c. To determine a calibration equation for calculating octane, 
aromatics, olefins, polyols, hydroxyl number, etc., for a near-infrared instrument, 
such as NIR Systems Model 6500, LT Industries Model 1200, Guided Wave 

1 1 Model 300 Series; Bran-Luebbe Infralyzer; or Infraprover; approximately 200 blends 

12 are used. Unfortunately, since gasolines are blends with components having 

13 different volatilities, it is very difficult, if not impossible, to store these samples to 

14 be run on an uncalibrated spectrophotometer at some later date. In other words, 

15 

the samples used to calibrate the reference spectrophotometer need not be 

17 available. Not using the same materials used to calibrate the reference 

18 spectrophotometer compounds the difficulty of calibrating an uncalibrated 

19 spectrophotometer because each member of the set of sample-s of the 

20 approximately 200 had to have its octane or knock values determined in 
accordance with ASTM D2699 and D2700, or the like. The procedure for 
determining knock values are very time consuming and require multiple 

24 determinations on the same sample with different engines before^ reliable 

25 knock-value results are obtained. 

26 Unknown prior to this invention was whether pure gasoline-like components, 

27 e.g. specific hydrocarbon molecules, that are readily reproducible could be 
selected for the set of sample-c. We have found in the case of near-infrared 
having wavelengths in the range 800 to 2500 nm (nanometers), that a set of 

31 sample-c comprising: benzene; methyl t-butylether; iso-octane; n-decane; 

32 1-octene; 2,3,4-trimethylpentane; n-heptane; cyclohexane; n-hexane; or toluene; 

33 when measured in both a reference instrument and an uncalibrated instrument to 

34 produce spectral data from each instrument can provide sufficient data or 
35 

information to calibrate the uncalibrated instrument. A calibrated instrument for 

36 

37 purposes of this specification is one that provides stable and reproducible spectral 

38 responses and has a calibration equation by which statistically significant values 

39 for physical or chemical properties of a sample can be inferred with an accuracy 
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equal to or exceeding the accuracy of a direct test method (e.g., ASTM) for 

5 determining such chemical or physical property. 

6 Broadly, to transfer a calibration equation for octane, any compounds that 

7 are near-infrared active, i.e. absorb near-infrared radiation, in 800 nm to 2,500 nm 

8 can be used. Examples of such compounds are: n-paraffins, iso-paraffins, 
9 

aromatics (such as dichlorobenzene or benzene), alcohols (such as methanol or 

n ethanol), olefins (such as 1-octene), halogenated alkanes (such as 

12 trichloroethane), ketones (such as methyl ethyl ketone), ethers (such as 

13 methyt-t-butyl ether), esters (such as amylacetate), and naphthenes (such as 

14 cyclohexane, and methyl cyclohexane). Throughout this specification and claims, 

15 

standard error of estimate or prediction are all at a 68% confidence limit, unless 

16 

17 expressly stated otherwise. Each is respectively referred to also as "SEE* and 

18 "SEP". 

19 The calibration equation for a particular instrument predicts a physical or 

20 chemical property based on spectral responses or spectral data. When measuring 

21 

octane values, the particular spectral responses or mathematical transforms 

22 

22 correspond to values of the spectral features that are selected for predicting 

2 4 octane values for a particular sample and are explained in detail in US 4,963,745. 

25 The equation is preferably linear, but non-linear equations may also be used. 

26 Spectral features are functionally related to spectral data, in that each 

27 spectral feature is uniquely defined by such spectral data. An example of such 
28 

functional relations are mathematical transforms such as derivatives of first, 

29 

3Q second, and higher order; fourier analysis yielding frequency components and their 

31 contributions; difference spectra; integration by segments; and the like. The only 

32 limit on choosing such functional relations or mathematic transforms is that any 

33 values derived by applying them to spectral data must be correlatable as 

34 independent variables to accurately determine the value for a dependent variable 
35 

in the calibration equation which corresponds, directly or indirectly, to the value for 

36 

37 a physical or chemical property. 

33 In discussing the invention in detail, the following definitions will be helpful. 

39 A "function" for purposes of this application is a correspondence that associates 

40 
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with each set of values for a set of independent variables, a unique value for a 

5 dependent variable. In accordance with this definition, a "calibration equation* of 

6 this invention is a function that establishes a correspondence between independent 

7 variables having values determined by spectra! features, and a dependent variable 
that corresponds to the value for a chemical or physical property. Any function 
that determines values for a spectral feature that are both uniquely associated with 

1 1 spectral data and can be used in a calibration equation is within the scope of the 

12 invention. Spectral data are uniquely defined by the results of measuring spectral 

1 3 responses of a material to electromagnetic radiation. Examples of electromagnetic 
radiation have been discussed in other portions of this specification. The prior art 
has established a variety of functions that can be used to derive values for spectral 
features from spectral data. So long as there exists a calibration equation which 

18 contains independent variables equal to the spectral features, defined by a 

19 function, there are no limitations as to the particular spectral features or functions 

20 for spectral data that determine such spectral features which can be employed with 

21 the invention. 
22 

"Calibration equation" for purposes of this specification and claims means 

24 an equation that predicts directly or indirectly; e.g., is functionally related to, the 

25 amount of a chemical or physical property (the dependent variable) based oh one 

26 or more spectral features (the independent variables) and perhaps some other 

27 parameters; e.g., temperature. 

A near-infrared regression model for octane measurements in gasoline 

29 

3Q which contain MTBE was published in a book of proceedings by the American 

31 Chemical Society, Division of Fuel Chemistry, in Volume 35, Number 1, Pages 266- 

32 275 which was written by Steven M. Maggard and also presented by him at the 

33 199th American Chemical Society National Meeting in Boston, Massachusetts on 

34 April 22-27, 1990. Specifically disclosed in this paper and presentation is the use 
of methyne t-butyl combination bands, including overtones, in models to predict 
octanes of MTBE-containing gasoline compositions. 

38 "Spectral features," for purposes of this application, means absorbance, 

39 transmrttance, reflectance (diffuse, specular, total, attenuated total, etc.), emission, 

40 
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fluorescence, phosphorescence, or transflectance, or mathematical functions 

4 

5 thereof such as derivatives, ratios, fourier transforms, and differences for ultraviolet, 

6 visible, near-infrared, mid-infrared, far-infrared, microwave, or radio frequency 

7 electromagnetic radiation. 

8 Throughout this text we have referred to the use of derivatives of spectral 
g 

absorbance data. Here we will show an example of what we mean when we use 

10 

1 1 the term derivative. The logic used to take a derivative of spectral absorbance 

12 data is based on the calculus of finite differences since the data is expressed at 

13 regularly spaced intervals and is discontinuous in between the spaced intervals. 

14 Spectral absorbance data such as gathered by a NIRSystems OL-6500 can 

15 be used. The OL-6500 has a spectral bandpass of about 9 nm but, the spectral 

16 

17 data is digitized such that recordings are displayed at 2 nm intervals. In general 

18 two things may be varied when the derivative about a point is calculated; the first 

19 of these is termed the segment and the second of these is termed the gap. The 

20 segment determines the number of points which are used to calculate the average 

21 number to the right and left of the point for which a derivative is being taken. The 
22 

gap on the other hand refers to the number of points about the midpoint 

24 (immediately adjacent to the point where the derivative is being taken) which are 

25 excluded from the calculation of the adjacent segments (although the total segment 

26 length is unaltered). If we refer to the average value of the segment (higher 

27 wavelengths) above the point at which the derivative is being taken as A and the 
average value of the segment below the point where the derivative is taken as B, 

29 

3q then according to the calculus of finite differences the first derivative of the point 

31 is A mints B. 

32 For example, when the segment is 10 in 2 nanometer increments/segments 

33 and the gap is zero, the first derivative of the spectral absorbance (d 1 A/dw) at 

34 1 170 nm is calculated as follows: 

35 
36 

37 d'A/dw^ = {(A l172 +A l174 +A l176 +A l178 +A 1180 )/5} - 

38 "f (^1 168 + 166 + ^1164 162 + ^116o)/^} 

39 If the segment is changed to six and the gap is changed to 8 the first derivative is: 

40 
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d A/dw 1170 = {(A, ,eo)/3> - {(A, 164 + A 1 18 2 + A, leo )/3} 



1 

2 
3 
4 
5 

6 In a general way, for data occurring at two nanometer intervals the first 

7 derivative of spectrai absorbance data at wavelength x is: 

8 
9 
10 
11 

12 -L (A 0t . tov/2) ^ 1) )/(segment/2) 

13 

14 where: 
15 



d'A/dw, = L (A & 

+ (o«p/2) *2q)/ (segment/2) 



x = wavelength in nanometers, e.g. 1170 and 
i = 1 to segment/2 



16 
17 
18 

19 In NlRSystems software package. NSAS (near infrared spectrai analysis 

20 software, NlRSystems, Inc.,), for reasons which are unknown, the above equation 
is strictly only true for segments of 2, 6, 10, 14, 18 etc., and for gaps of 0, 4, 8, 12, 
16, etc. A segment of 0 is treated as a segment of 2. Lfeewise a segment of 4 is 

24 treated as 6, etc. Similarly a gap of 2 is treated as a gap of 0, and a gap of 6 is 

25 treated as a gap of 4, etc. 

26 Higher order derivatives may be simplified to expressions involving only first 

27 derivative terms and solved by the form of the previous equation. In general when 
the derivative being sought is of the h-th order: 



30 h 

31 d A/dw x = d^A/dWix+Ksegmert+gttf/ty+i) 



33 " d WdW^Mgnwit+BapJ/aM). 

34 



where the gap is equal to gap of the h-th ordered derivative and the segment used 
to calculate the (h-1) -th ordered derivative is the same as that being sought for 
33 the h-th ordered derivative. 
39 
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4 For example, the 4th derivative of 11 70 nm at a segment of 10 and a gap 

5 of 4 is: 
6 

7 d 4 A/dw 1170 = dtydw,^ - d 3 A/dw 1162 
8 

^ = (d 2 A/dw 11fl6 - dWcKiroM^A/dw^ - d 2 A/dw 1164 ) 
11 

12 = {(dWdWnw - d 1 A/c»w 1178 )-(d 1 A/«lw,i78 - d 1 A/dw 1162 )} 

13 

14 -{(d'A/dw,^ - d 1 A/dw 1162 )-(d 1 A/dw 1162 - d'A/dw,,^} 
15 

17 = {{[((A 1198 + A 12M +...+A 1206 ) /5H(A 1190 +A 1188 + +A 1182 )/5)J 
18 

1 9 - t((A 1 1 82 + A 1 1 84 + • " + A 1 1 90 )/5) - ((A 1 1 74 + A 1 1 72 + - + A 1 1 66 ,/5)1> 
20 

21 •<I« A 1182 +A 1184 + - +A 1190 )/5W(A 1174* A 1172 + " 4A 1166 )/5)J 
22 

23 

24 l(tA 1 166* A 1 168 + " +A 1 174>^"H^1 158* A 1 156* "* A 1 160 )/5)] » 

25 

26 • {{l( ^1182 tA 1184 + " +A 1190 ,/5H ^1174 +A 1172 + " +A 1166 ) ^ 
27 



28 
29 
30 



< «< A 1 1 66 + A 1 1 68 * ■ * A 1 1 74 ,/5HtA 1 1 58 * A 1 1 56 * " +A 1 1 SO*' 5 " 

29 

<W<\ 166 +A 1 168 + " +A 1 174 ,/SH(A 1158 +A 1 156* - +A 1 150 )/S,J 



31 
32 

33 f ^1150 +A 1152 + " +A 1158 >/5H ^142 +A 1140 + " +A 1134 )/w > 
34 

35 

where all of the lower ordered derivatives are also taken using a segment of 10 

36 

37 and a gap of 4. Note that it is necessary to correct the segment and gap to a 
33 valid value before applying one of the generalized formulas (e.g. the fourth 

39 
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3 

4 derivative at a segment of 8 and a gap of 6 is corrected to a segment of 10 and 

5 a gap of 4 before applying the formula). 

6 Note segments of 2 nanometers were used, but segments of smaller or 
larger numbers or fractions of nanometers can be used. Tne choice preferably 
depends on the degree of spectral resolution available to the spectrometer. 

The Examples of differentials, finite differentials, discussed are but one of 
, , many of the types of finite differentials that may be taken. The theory behind all 

12 of the methods that are available are substantially the same; however, the specific 

13 details involved in the calculation may vary somewhat. It Is also to be noted that 
although finite differentials may be discussed, it is possible that differentials having 
differences that approach zero, rather than some finite number are also possible. 

17 For example, a polynomial power series like a Taylor Series can be established 

18 based upon a series of points taken from intensity versus wavelength curves or 

19 spectral data. 
This invention provides a means for determining a calibration equation for 

a primary instrument that takes advantage of a substantial amount of work 
previously expended to obtain a calibration equation for a reference 

24 spectrophotometer. The primary instrument may, in fact, be the reference 

25 spectrophotometer that has gone out of calibration for reasons that have already 

26 been discussed. ; ; 
To make the discussion easier to understand, the following conventions will 

be followed. The independent variables appropriate to a calibration equation for 
a primary instrument will be identified by H" or "-1", e.g., independent variable-i 

31 or independent variable-1 corresponds to the one or more members of a set of 

32 independent variables that are contained within the calibration equation, calibration 

33 equation-i, or calibration equation-1, respectively, for the primary instrument, 
referred to throughout as instrument-i or instrument- 1 . in most examples, there are 
three wavelengths used, each wavelength is referred to as wavelength-j, from 
which three independent variables-i derive values functionally from spectral 

3 8 responses of a sample, c, taken from a set of samples referred to as sample-c. 

39 Corresponding to a particular set of independent variables-i, there is a dependent 
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variable which will be referred to as dependent variable-i. The value of dependent 

4 

5 variable-i will usually be equal to a predicted value of property-p of sample, c. 

6 Dependent variable-i corresponds to that dependent variable found in calibration 

7 equation-!. The constants in calibration equation-i are referred to as constants-i. 

8 There are two different types of constants-i. One corresponds to a slope constant- 

g 

i, and the other to a bias constant-i. The slope constants-i are coefficients of the 

10 

1 1 independent variables-i. Clearly, there is a slope constant-i for each independent 

12 variable-i which is in turn dependent upon the number of wavelengths-j that are 

13 used in the calibration equation-i. There is only one bias constant-i which is not 

14 a coefficient of any independent variable-i. 

15 When one determines a calibration equation appropriate to instrument-k 

16 

• (also referred to as instrument-2) by recourse to the calibration equation for 

18 instrument-i, one essentially determines a new set of calibration equation constants 

19 that correspond to constants-k. New constants are found in at least one and 

20 preferably two steps, assuming that one already has a calibration equation-i 

21 already determined for instrument-i. 

22 

To successfully transfer or modify calibration equation-i into calibration 

23 

24 equation-k, the SEE or SEP of values for property-p based upon spectral data of 

25 instrument-k or instrument-2 by means of calibration equation-k or calibration 

26 equation-2 should most preferably be substantially as good as that achieved by 

27 instrument-i and calibration equation-i or alternatively be substantially as good as 

28 that which would be achieved had instrument-k been calibrated the same way as 

29. 

instrument-i, but in most cases be at least equal or better than the SEE for the 

30 

31 primary test method for directly measuring values for property-p used to calibrate 

32 instrument-i. The benefit made available by this invention is that the number of 

33 samples that must be run on either instrument-i or instrument-k, or both, is 

34 significantly less than that number of samples required to obtain a reliable 

: * ' 35 

calibration equation-i for instrument-i or to fully calibrate instrument-k in the same 

37 98 instrument-i. In other words, whatever the value, E(k), in the standard error 

38 °* prediction for the value for a property-p based upon instrument-k and calibration 

39 equation-k will not be substantially greater than the value for the standard error of 
40 
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4 prediction resulting from determining values for a property-p by the direct or 

5 primary test method initially used to calibrate instrument-i. 

6 One way, but not the only way, to transfer calibration equation-i to 
instrument-k, is to change the constants of calibration equation-i. Tnese changed 
constants become constants-k. Calibration equation-i is preferably transformed 

iq int0 calibration equation-k, in two steps. Step Two (below) is optional, but 
n preferable depending upon whether the correction indicated by step two is 

12 significant. By significant, is meant the inferred variance in predicted values is not 

13 more than about 10% greater than SEE or SEP, whichever is the greater, achieved 
by calibration equation-k without Step Two. 

Step One: Determining calibration equation-k (calibration equation-2) by 

17 means of calibration equation-i (calibration equation- 1) involves the following 

18 sub-steps: 

1 9 (a) Determine calibration equation-i by measuring dependent variables-i 
directly by a primary test method with a standard error of estimate 
of E, and correlating said dependent variable-i to a set of 

~ independent variables-i. The set of independent variables-i may 

24 consist of one or more members. ' 

25 (b) Obtain respectively with instrument-i and instrument-k spectral data-i 

26 and spectral data-k, respectively, for each member of sampte-c, 
which are to be used in step (c) (comprising at least one member, 
c); functionally determine value-k and value-i, respectively, (equal to 
each spectral feature-k and -i, respectively) for each member of 

3! sample-c corresponding to each independent variable-k and -i, 

32 respectively. 

33 (c) Determine (in accordance with calibration equation-i) a value for said 
dependent variable-i for each member of sample-c. The value 
determined for dependent variable-i or a value determined by some 

37 function-f using that value determined for dependent variable-i is 

38 inserted into calibration equation-k for the value of dependent 

39 variable-k in calibration equation-k. Also insert into calibration 
40 
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equation-k the values for spectral features-k determined for each 
member of sample-c. (Recall that spectral features-k are functionally 

6 determined by spectral data-k that is determined from measurements 

7 made by instrumerrt-k of each member of sample-c. The only 

8 unknowns in calibration equation-k will therefore correspond to the 
constants-k.) The values for constants-k that are to be substituted 
into calibration equation-k produce a resulting calibration equation-k, 

12 wherein the sum of all absolute differences between each estimated 

13 value of property-p for sample, c, corresponding to dependent 

14 variable-k for sample, c, and that estimated value of property-p for 
sample, c, corresponding to dependent variable-i for sample, c, is 
minimized so that values for property-p for sample-u that have values 

18 for property-p within the range of interest are inferred based upon 

19 the resulting calibration equation-k, with an SEP of not substantially 

20 greater than the SEE or SEP of E for determining dependent 

21 variable-u directly by some primary test method. Once one has 

defined the constants-k, that permits one to infer the value of 
dependent variables-k for members of sample-u, having values for 

25 property-p, within the range of interest, within the SEE or SEP of E 

26 of the direct or primary test method, then one has obtained a 

27 calibration equation-k in light of the former calibration equation-i. An 

28 absolute difference is the value for a difference betwe n two 

29 

numbers that results from first squaring and then taking the positive 

30 

31 square root. (Recall that instrumerrt-k, the primary instrument, may 

32 in fact be the reference instrument that has gone out of calibration, 

33 or it may in fact be a totally different instrument). 

34 The value for dependent variable-i that is to be inserted into 
calibration equation-k for dependent variable-k can be functionally related 
to dependent variable-i by a function-f such as for example a logarithm [log) 

33 of any base (e.g., natural or base 10) such as Y e k = nelog Y,,' + m or as 

39 a linear function, Y e k = n Y c '+m, where n and m are constants, or a linear 

40 
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1 
2 
3 

4 combination of these; or in general, any function such as a polynomial or 

5 geometric series; provided that over the range of values of interest for Y,", 

6 the function is substantially linear. 

7 The range of values of interest for Y, fe is that range of values for a 
chemical or physical property that calibration equation-k in conjunction with 



10 instrument-k have been calibrated to predict with an SEE or SEP of less 

j 1 than that SEE or SEP of the primary test method. 

12 Clearly, if Y k 0 = neV e + m and n = 1 and m=0, then Y k « Y 1 . 

1 3 However, for any other values of n and m, the corresponding constants will 
be related to those otherwise determined for n=l and m=0 as follows: 

A k - nV Bj k - nB )0 k 

17 where A," ■ value of constant for n=l and m=o and 

18 B Jc k « value of constant for n«1 and m«0 at wavelength-j 

19 ,n other words, whatever values for n and m that are used, the dependent 
variable-k will be functionally related to dependent variable-i. To determine the 
value for the corresponding physical or chemical property directly measured, one 
merely transforms inversely from Y c k to Y n k where for any m and n, except n=0, 

24 V = (Y c k -m)/n; " 

25 lnstrument-k and its spectral data are useable in this invention, if 

26 instrument-k could be calibrated by an identical method to that which could be 

27 used to successfully calibrate instrument-i. 

28 

2g Preferably, calibration equation-i and calibration equation-k are linear. In the 

30 generalized method discussed above for this invention, calibration equation-i and - 
3.1 k n ©«d not either be of the same form or both be linear. By the same form is 

32 meant having the same functional dependence on the same number of 

33 independent variables with corresponding constants determining a dependent 

34 variable. 

Preferably, calibration equations for instrument-i and instrument-k have the 
following linear form: 

38 Y.' - A'+I,{B, , X j8 , >andY. k = A k +I J {B j k V} 

39 where: 

40 
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4 Y, 1 Is the dependent variable-i in calibration equation-! that preferably 

5 assumes a value substantially equal to property-p of any member, s, of said 

6 sample-s and Y a k is the dependent variable k in calibration equation-k which 

7 is functionally related to Y, 1 by a function-f; 

8 A k and A* are each intercept or bias constants-i and -k in calibration 

9 

10 equations-i and -k, respectively; 

1 1 Xj, 1 and Xj, k are each one of said at least one independent variable-i 

12 and -k in calibration equations-i and-k, respectively, having a value equal to 

13 at least one said spectral feature-i and -k, respectively of a member, s, at 

14 wavelength-], and; 

15 

Bj and B, are slope constants-i and -k, respectively, corresponding 

16 

17 to coefficients for each said at least one independent variable-i and -k 

18 corresponding to spectral features-i and -k at wavelength-J, respectively. 

19 There may be virtually an infinite variety of ways of substituting values for 

20 constants-k that yields a calibration equation-k with an SEP not substantially 

21 greater than the SEE or SEP necessarily present in the direct method used at 
22 

some point to obtain directly or indirectly calibration equation-!. By "substantially 

24 greater is meant 30% greater, but preferably, 20% greater, and still more 

25 preferably, 10% greater. Directly determined values of property-p must at least 

26 once be obtained by some primary or direct test method, in order to calibrate an 

27 instrument by determining Its calibration equation for predicting values for a 

28 

physical or chemical property, property-p, for an array of materials, based upon 

29 

spectral data for such materials* However, once one has calibrated a first 

30 

31 instrument one can use such instrument to calibrate other instruments without 

32 necessarily remeasuring property-p by the primary test method, Instrument-i has 

33 a directly-determined calibration equation-! if some or all of the constants-i therein 

34 are determined by values directly determined by the primary test method. 

35 

However, where only values used to determine constants-k arise from values 

36 

37 predicted by a previously calibrated instrument, without values from the primary 

3 8 test method for property-p, then instrument-k is indirectly calibrated. In all cases, 

39 there will be at least one instrument whose calibration equation is determined 

40 
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4 primarily based upon values for its dependent variable determined by recourse to 

5 a primary test method. It is the SEE or SEP of that primary test method which is 

6 being referred to above. 
However, we have found that the following criteria are particularly useful in 

determining constants-k assuming that calibration equations-k and -i have a linear 
form as defined above. One criterion requires that constants-k correspondhg 
respectively to B" and A K minimize to substantially the maximum degree possible 

12 the following expression: £ e (Y e k -A k -I) J B | k X Jc k ) 2 . 

1 3 Alternatively, values of B, k and A k can be used which substantially satisfy the 
following two equations simultaneously: B, k = BjEJ^ X^/^q^ 2 and A k = 2 c (Y e k - 
2 rWV n ; where n - number of members of sample-c. 

17 Finally, we have also found a particularly useful criteria for determining 

18 constants-k corresponding to Bj k and A k is one which involves substantia^ 

19 satisfying the following two equations simultaneously: Bj k = B j 'E c X ie l /E c (X |c k ) and 

20 A k = E e (Y c k -E J B j k X je k ) / n; where n = number of members of sample-c. 
Step Two involves a slope and bias correction. The requirements for 

samples appropriate to a slope and bias correction are discussed in the Examples. 

24 Essentially, a slope and bias correction arises from determining values for "m" and 

25 "b" which satisfy the following expression: 

26 m = KIcY 0 a £Y e Vn-£(V YjA/KEoY." ZeY 0 k }/n-Z e Y e k Y k ] 

27 b - {E e Y e *-mI 0 Y c k }/n 

28 where: 
29 

30 Y e is a directly or indirectly determined value for property-p of a 

3T particular member, c, of sample-p that is determined with an inferred 

32 standard error of estimate no greater than E of the primary test method; 

33 and wherein said value for property-p is within the range of interest; 

n is the number of members of sample-p; and 
v is the standard error of estimate from using calibration equation-k 

37 values to determine dependent variable-k, for sample member, c, wherein 

38 v is not significantly greater than E. 
39 

40 
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Substitute said values for m and b into the following expression and 

4 

5 an improved calibration equation results for the dependent variable-k for 

6 instrument-k, which is as follows: 

7 Y e " = [mA^b]+S[mBfV]. 

8 Step Two may in some instances be omitted because the correction 

9 suggested or required by applying a slope and bias correction does not result in 

1 1 a significant change in constarrts-k for calibration equation-k. But it is preferable 

1 2 to carry out the steps involved in a slope and bias correction to determine whether 

13 a significant correction is involved or not We have found in most instances, a 

14 slope and bias correction is required. 

15 What is surprising about Step Two, the slope and bias correction, is that 

16 

17 even a significant variation in values of constants-k required by or due to the 

18 correction does not mean that an error has occurred in carrying out Step One. 

19 We have discovered that the error, if any, in Step One is substantially corrected by 

20 a slope and bias correction and vice versa. (See Example 1). Further, we have 

21 found that the number of samples required to carry out a slope and bias 

22 

correction, even when added to those required to carry out Step One, are together 

23 

nA far fewer in number than was required to determine calibration equation-i for 

24 

25 instrument-i. This is totally unexpected and could not have been anticipated. In 

26 fact, when a significant correction in constants is indicated by a slope correction, 

27 the art (see manual for Near Infrared Spectral Analysis Software published by NIR 

28 Systems, Inc. Silver Spring, Maryland, 1989) teaches that one should recalibrate 

29 

instrument-k by the same method used to initially calibrate instrument-i. 

30 

31 Closer examination of what is happening to the constants-k in calibration 

32 equation-k as one applies Steps One and Two is interesting. In Step One, each 

33 constant varies substantially independently in order to satisfy the set of numerical 

34 criteria discussed above. In Step Two, slope constants all vary together by an 

35 amount determined by "m". Intercept constant-k is impacted in two ways, both 
36 

by *m" and an incremental amount *b". When applied together, a calibration 

37 

33 equation-k for instrument-k preferably provides an SEP for dependent variable-k 

39 that is almost as good as, and often much better than that achieved by calibration 

40 
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3 

4 equation-i using data measured by Instrument-i or at least an SEP that is equal to 

5 or better than E, the value of SEP for the direct test method for determining values 

6 used to calibrate instrument-i. This result shown in the examples, is totally 

7 unexpected, surprising and valuable. 
Also contemplated by this invention is the indirect method for determining 

a property-p for a sample-u using spectral data-k measured by an instrument-k, 
1 , capable of determining such spectral data-k for sample-u. Though the art teaches 

12 many ways for identifying a calibration equation for ah instrument, the art has not 

13 provided up till now an easy method for transforming a calibration equation-i 
appropriate for instrument-i to a calibration equation-k appropriate for instrument-k. 

A unique form for calibration equation-k or -i involves the use of adjacent 

17 order differentials as spectral features. For example, the calibration equation for 

18 instrument-k and instrument-i each have the form: 

19 Y^A'+T^V+^V 

20 Y^A'+^Bj'^'+T^q'Xj,," 

Here Y e ', and Y e k are calculated or predicted values for dependent variables 
based upon the calibration equations for instruments-i and -k, respectively, both 

24 for each member, c, from sample-c. 

25 Xjc 1 ' and X|c k are second derivatives of observed absorbances for 

26 instruments i and k, respectively, both for each member, c, from sample-c at 

27 wavelength-j. 
Xje 1 , and Xjc* are the third derivative for observed absorbances for 

instruments-i and -k, respectively, both for sample-c at wavelength-j. 

31 We have discovered that with calibration equations in the form of adjacent 

32 orders of derivatives of spectral data, that surprisingly a slope and bias correction 

33 brings about a variation in calibration equation constants-k that results in inferred 

34 values for property-p of a sample, c, based upon spectral data measured for 

35 

sample-c by instrument-k that is in excellent agreement with the results found from 

36 

37 instrument-i and calibration equation-i. 

38 
39 
40 
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3 EXAMPLE 1 
4 

5 A calibration equation for an uncalibrated instrument is obtained based upon 

6 the calibration equation established for a reference instrument. Solutions for the 

7 following equations determine: B, k, s and A k in the inear calibration equation for 

8 instrument-k, the uncalibrated instrument, having the form: Y e k ■ A k +£B, k X Jc k ; 
which relies on calibration equation: Y e ' = A' + ZjB/Xj,'; for a calibrated or reference 

j j instrument, instrument-i. 

12 Three different mathematical criterions are used Independently. Since only 

13 three wavelengths are used, J = 1 to 3; and since 25 = total number of members 

14 in sample-c, c = 1 to 25. Independent variables, X^, (where "k" identifies 
instrument-k) have values for spectral features derived from the spectral response 
of instrument-k at wavelength-j. Each wavelength-j is within the wavelength range 

18 of from 800 to 2500 nm. Here the spectral features are the second derivative of 

19 the absorbance determined at three wavelengths-J: 1196, 1220, and 1286 nm. 
20 

21 Mathematical Criterion I 

Io( Y o , -A k -Z,B l % k ) = 0; Assume Y e ' = Y 0 k 

IcYc V-A'E^^^EPC,^ = 0; 

25 2UY e V-A k EeX^B>2^^ = 0; and 

26 &Y e V-A>I^ k -B>Z c X^^ * 0 

27 Instrument-i is caibrated using multiple linear regression at the second 

28 derivative of absorbance at 1196, 1220, and 1236 nm modeling the dependent 

29 

variable (R+M)/2 octane using 72 gasoline samples encompassing the octane 

30 

31 range of 86-93 in an approximate box car distribution. Instrument-i is a NIR 

32 Systems 6500 operating in transflectance mode using fiber optics at a path length 

33 of 18mm. 

34 From the regression equation for instrument-i, one obtains values of 25 

35 calibration transfer standards and subsequently assigns these as values of the 

36 

3? dependent variable (R+MJ/2 octane for the 25 calibration transfer standards. The 
38 standards are the same ones seen in Example 2. 

39 
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3 

4 The 25 calibration transfer standards and 10 samples of known octane are 

5 run on instrument-k. The results are summarized in Table 1 . 

6 Alternatively, we have also determined that it is possible to indirectly derive 
calibration equation-k from the caiibrabon transferred equation-k, where a least 
squared error reducing step was used in place of the slope and bias adjustment 
When the errors were modeled as the dependent variable-k using the spectral 

n data-k, again as the independent variable, it was possible to sum the resulting 

12 constants with the constants from the calibration transferred equation-k and to 

13 indirectly derive the least squared solution (mathematical criterion 1). 



Mathematical Criterion II 
B, k = IB^Vl/tlo^) 2 ] 

A k = k(Y e k -LBX k )/n A§iumt Y> , m y k 



NOTE: c = 1-25; since only three wavelengths 1196, 1220, and 1236 are used; 
and n = 25; since number of members In sampte-c Is 25. 

22 

23 lnstrument-i is calibrated using multiple linear regression at the second 

24 derivative of absorbance at 1196, 1220, and 1236 nm modeling the dependent 

25 variable (R+M)/2 octane using 72 gasoline samples encompassing the octane 

26 range of 86-93 in an approximate box car distribution. Instrument-i is a NIR 

27 Systems 6500 operating in transflectance mode using fiber optics at a path length 

28 of 18mm. 

29 ^ _ 

30 From the regression equation for instrument-i, one obtains values of 25 

31 calibration transfer standards and subsequently assigns these as values of the 

32 dependent variable (R + M)/2 octane for the 25 calibration transfer standards. The 

33 standards are the same ones seen in Example 2. 
The 25 calibration transfer standards are run on instrument-k. Regression 

constants are calculated as follows: 
B i k ■ B i2» k l B i22o Zc X 12a0ie ' X 122fte k ] / £p (X^^*) 2 
38 Ak s le IY e ' - (B^^O - (B^O - B^X,^*)]/^ 
39 
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3 

The above constants in the calibration equation predict the value of 10 

4 

5 sampLs of known octane. There is a severe bias of 5.314 (R+M)/2 octane 

6 numbers. A computer program establishes that a stops correction of 0.56741 is 

7 needed and an intercept correction of 35.646. These adjustments yield the final 

8 slope and bias adjusted calibration equation for instrument-k and show a standard 
9 

iq error of estimate of 0.184 and bias of 0.0265 (R+M)/2 octane numbers on the 

11 remaining 62 samples. 

12 

13 Mathematical Criterion lil 

'< = B f 'T^'/IeV 

A k = Ee^-IrBrX")/" Assume Y c ' = Y c k 

NOTE: r = 1-3; since only three wavelengths are used, n = 25; since number of 

18 members in sample-c is 25. Just as with Criterion II, the above Criterion III are 

19 used. The results of all methods are given in Table 1 . 

20 To demonstrate the relative accuracy to be expected from an octane 

^ calibration equation transferred from one instrument, instrument-i to another, 
22 

instrument-k, the following is done. 

23 

24 Both instrument-i and instrument-k are each independently calibrated using 

25 the same 72 samples. The calibration establishes a linear calibration equation that 

26 can be used to predict octanes. A standard error of prediction for each of these 

27 equations Is established for ten knowns. 

28 The transfer calibration equations for instrument-k is derived from calibration 

29 

equation-i and one of the following criterion: I, II, or III. 

31 Each transferred equation is slope and bias corrected with ten samples 

32 yielding a slope and bias corrected calibration equation. 

33 The results from the above steps are shown In Table 1 . What is particularly 

34 significant is that the SEP dramatically changes as a result of the slope and bias 
correction. However, the final SEP, after performing a slope and bias correction, 
in the case of Criterion I, yields an SEP identical to that of the original calibration 

38 equation-i. 
39 
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A good estimate of the best SEP one might expect in this instance would 
be an SEP of 0.146 which is the SEP for instrument-k when fully calibrated using 
the 72 samples. 

TABLE 1 

Calibration Transfer Without Slope & Bias Correction 



r 


Trar 

Equj 
Stani 


isfer Calibration 
rtion-k Using 26 
dards of Table 1 


Original 
Calibration 
Eauatlon-I 
(72 Samples) 


inavperraeni 
Calibration 
of 

Ifififn immil fc 
%m u ui irfjf i%*fk 

(72 
Samples) 


Criterio 
n I 


Criterion 
II 


Criterion 
III 






A* 


74.501 


76.041 


75.560 


76.263 


70.918 


B, h 

(1220nm) 


134.620 


206.02 


211.90 


124.60 


131.21 


B k 

(1 2 196nm) 


21.139 


26.939 


26.940 


24.407 


13.48 


Ba k 

(1236nm) 


69.345 


70.681 


70.689 


65.546 


65.251 


SEP 


0.366 


5.383 


5.523 


0.156 


C ai46 


Bias 


-0.288 


5.195 


5.314 


0 


0 
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Calibratlon Transfer With Slope & Bias Correction Using 10 Sample* 





Predicted Results of 62 
Unknowns Using Transfer 
Calibration Equation* 


Independent 
Calibration 

01 

Instruments 
(72 

Samples) 




Criterion 
I 


Criterion 
II 


Criterion 
III 






76.013 


78.792 


78.753 


70.918 


B," (1220nm) 


123.20 


116.90 


116.86 


131.21 


I B 2 * (1196nm) 


19.345 


15.285 


14.857 


13.48 


B 3 k (1236nm) 


63.460 


40.105 


38.985 


65.251 


SEP 


0.156 


0.184 


0.188 


0.146 


Slope Adjustment 
(m) 


0.91514 


.56741 


0.55150 


1 


Bias Adjustment 
<P) 


7.8341 


35.646 


37.082 


0 


Bias 


0.0191 


0.0265 


0.0373 


0 



EXAMPLE 2 

In this example, Mathematical Criterion I of Example 1 (multiple linear 
regression) is used and the form of the calibration equation is Rnear, but spectral 
features equal to second and third adjacent derivatives of spectral data 
corresponding to each selected frequency are used for the independent variables. 

Y c k =A k + I j B j % k + I je C j k X ie ,k 

Y.'-A'+SB" Xj.'+^C/ V 

Here Y c ', and Y c k are calculated or predicted values for dependent variables 
based upon the calibration equations for instruments-i and -k, respectively, both 
for each member, c, from sample-c having 8 members. 

Xje', and X^" are second derivatives Of observed absorbances for 
instruments i and k, respectively, both for each member, c, from sample-c at 
wavelength-j, where j = 1 to 3. 
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3 * 3V 

4 X^", and Xp* are the third derivative for observed absorbances for 

5 instruments-i and -k, respectively, both for sample-c at wavelength-j, where 

6 j = 1 to 3. 

7 The three wavelengths are the same as used in Exampie 1. 
Seventy-two gasoline samples are analyzed on instrument-i to obtain their 

absorbance spectra. In addition, 25 samples of calibration transfer standards are 
n also run. The second and third derivatives of the spectra are then made using 

12 segment lengths of 20 and gaps of 0 on NSAS (Near Infrared Spectral Analysis 

1 3 Software, available from NIR Systems, Inc.). The compositions of the 25 calibration 
transfer standards are shown in Table 7, and are specifically chosen because they 
encompassed the majority of variance of the 72 gasoline samples at the 

17 wavelengths used in the calibration. The 72 gasoline samples have octane in the 

18 range of 86.3 - 93.2 (R+M)/2 octane numbers and are distributed among this 

19 range in an approximate box car distribution. Each sample is analyzed four times 
by four different knock engine laboratories (APAL, No. 2 Control Lab, Canton 
Control Lab, and St Paul Control Lab) by ASTM methods 2699 and 2700 to 
provide accurate estimates of (R+M)/2 octane. 

24 The second and third derivative files of the 72 gasoline samples are sent to 

25 a mainframe computer, and a multiple inear regression in accordance with 

26 Mathematical Criterion I is performed by conventional computer programs using 
SAS statistics (available from SAS Institute, Inc., Cary, North Carolina) on the 
second derivative of absorbance at 1196, 1220, and 1236 nm and the third 
derivative of absorbance at 1196, 1220, and 1236 nm modeling the dependent 

31 variable of (R+ M)/2 octane. Analysis of the data in this way provides a multiple 

32 linear regression model of the form: 

33 

Y e '=[(R+M)/2J octane = A^B^+^'p^VB.'^+C^M+Q.'^+C,' PO 

34 
35 

36 The values of the B, terms are shown in Table 2. The multiple coefficient of 

37 determination (R 2 ), multiple coefficient of correlation (R), F-value and standard error 

38 of estimate are shown in Table 2. It can be concluded from Table 2 that the 

39 calibration equation is extremely accurate. 

40 
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TABLE2 

Regression Constants Obtained f r the 
Calibration Equatlon-l 
of the 72 Gasoline Samples 



Variable 


Wavelength 


Derivative Used 


Value 1 


A' 


Constant 




71.18951 I 


b; 


1220 nm 


2nd 


46.25463 I 


B 2 ' 


1196 rtm 


2nd 


-22.35057 | 




1236 nm 


2nd 


64.61244 I 


c,' 


1220 nm 


3rd 


-28.34514 I 


c 2 ' 


1196 nm 


3rd 


-12.83303 I 


c 3 ' 


1236 nm 


3rd 


-66.41106 I 



Statistical Results For the Calibration Equation^ 



| Multiple 


Multiple 




Standard 


R Coefficient of 


Coefficient of 




Error of 


I Determination 


Correlation 


F-Value 


Estimate 


1 0.9962 


0.9980 


2805 


0.145 



Next the values of the Calibration Transfer Standards are predicted from the 
calibration equation of instrument-i and the predicted values are used as the 
dependent variable ((R + M/2) octane) and the absorbance data for the Calibration 
Transfer Standards on instaiment-k are converted to second derivative and third 
derivative and are used for the values of the independent variables-k to determine 
relevant constants-k using the classical least squares method or Mathematical 
Criterion I. The results are given in Table 3. 
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TABLE 3 

Regression Constants Obtained on Instrument-k For 
the 25 Calibration Transfer Standards 



Variable 


Wavelength 


Derivative Used 


Value 


A" 


Intercept 




55.55131 


B," 


1220 


2nd 


178.49118 




1196 


2nd 


-47.22965 




1236 


2nd 


130.60840 




1220 


3rd 


-71.70553 


c." 


1196 


3rd 


-50.71056 


of 


1236 


3rd 


18.43823 I 



Table 4 shows the statistical data for calibration equatiook. 

TABLE 4 

Statistical Data for the Calibration Transfer 
to instrument-K Using the Octanes 
Predicted from Instrument-! 



| Multiple 


Multiple 




Standard ^ 


1 Coefficient of 


Coefficient of 




Error of 1 


I Determination 


Correlation 


F-Value 


Estimate 1 


1 0.9988 


0.9994 


2544 


0.124 | 



Seventy-two gasoline samples were then run on instrument-k and second 
and third derivative files are obtained as before, the new calibration equation from 
Table 3 is used to predict the (R+M)/2 octane of the 72 gasoline samples [note 
that these samples are identical to the original calibration set run on instrument-i] 
using the same octane as is used in the original calibration of instrument-i. These 
calculations are done on a PC using a BASIC program. The calculated octanes 
show a bias of -0.431 and have a standard error of prediction of 0.586. These 
errors are quite high considering that the original standard error of estimate on 
instrument-i is only 0.145. Also, it could be seen from the predicted octanes that 
the error for low octane predictions ((R+M)/2 approximately equal to 86) is quite 
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hlgh 0.e., about -0.8) when compared to the higher octane ((R+M)/2 
approximately equal to 93) predictions (e.g., much closer to zero, about 0.1). This 
type of behavior is similar to heteroscedasticity, in that calculated predictions tend 
to be erroneous by an amount that is a function of estimated octane. For 
example, lower octanes yielded high errors and higher octanes yielded lower 
errors. When the above behavior is observed, a slope and bias adjustment is 
warranted providing that a plot of (Y actual - Y predicted) versus octane 
approximates a straight line. This ordinarily is the case providing that the original 
calibration equation is linear or a non-iriear equation transformed so as to obtain 
a linear relationship. From examining the data, it Is clear that a slope adjustment 
is warranted since the relationship approximates a straight line. A PC 
BASIC program Is used to accomplish this and shows that a slope correction for 
"m" of 0.86406, and a bias correction for "b" of 12.5590 is needed. 

Y(Predicted) = Y c k = [0.86406(A k +L-i 72 B i k X jc k +L = , 72 C, k X^)]* 12.5590 

Table 5 shows the final regression constants for instaiment-k using the 
calibration transfer followed by the slope and bias correction. 

This equation is found to have a standard error of estimate of 0.190 on 
instrument-k compared to instrument-i's standard error of 0.145 on the 72 gasoline 
samples. 

TABLE 5 

Final Regression Constants for the Slope 
and Bias Corrected Transferred 
Calibration for Instrument-! 



| Variable 


Wavelength 


Derivative Used 


Value 


I A< 


Intercept 




60.55866 


B,' 


1220 


2nd 


154.22709 


B * 


1196 


2nd 


-40.80925 


B '' 


1236 


2nd 


112.85349 




1220 


3rd 


•61.95788 


c; 


1196 


3rd 


-43.81697 


1 c,' 


1236 


3rd 


15.93174 



r- 

► 94/08225 PCT/US92/08547 
1 

2 -35- 
3 

4 EXAMPLE 3 

5 This example demonstrates the slope and bias correction referred to earlier 

6 m this specification (Analytlca Chemica Acta, 192(1 987), 197). 

7 The various terms are as previously defined throughout this specification, 
the Y-term is a dependent variable calculated or predicted based upon 

1 Q independent variables equal to one or more spectral features identified as X-terms. 

1 1 Depending upon the size or magnitude of the slope and bias corrections 

12 corresponding to m and b; where m is the slope correction and b is the bias 

13 correction, one may decide not to modify the calibration equation for the 

14 uncalibrated instrument. Y k is the slope-corrected and bias-corrected value 

15 

predicted for a physical property as shown hereinafter. 

16 

17 In the specific example Y e \ the predicted octane value from calibration 

18 equation-k, before the slope and bias correction, and Y e \ the sample measured 

1 9 octane (or octane predicted by a calibrated NIR instrument) , are used to determine 

20 the slope and bias corrections. In the sample set that is used to do a slope and & 

21 bias correction, it is necessary that the set * 

22 

23 1- Be representative of the range of interest for a particular property, for 

24 example, octane. The sample then should have actual octane values 

25 (measured in accordance with ASTM D2699 and D2700) approximately in 

26 the range of octane values of interest. 

27 2. Have the octane measured for each sample, and that this actual octane 

28 

value be then correlated (in effect plotted) versus the calculated octane 

29 . , 

3Q value for each sample determined by using the calibration equation to which 

31 the slope and bias correction is to be made. 

32 3- The number, n, of samples used for making the slope and bias correction 

33 must be 2 to define both m and b, but in the examples, 10 is usually used. T 

34 Y e k = A k +£ e B | k X je k The calibration equation for instrument-k 

35 

Y e * is the directly determined, or actually measured property; here ASTM 

36 - 

37 measured octane. 

38 
39 
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4. To find m and b, solve for m and b which satisfy the following mathematical 
criteria: ^(YcVJ^^cX"^^ s 0 

6 LW-mY^b) = 0 

7 The results are: 

8 m = [(LeY^ EYoVn-IoCYo' Y c k )]/[{L»Y c k J^Y p k }/n^f p k Y c k ] 

9 b = IMIcY^-m^} 
Y c k - [mA'^+bl + IjfmB,^] 



10 



12 

13 EXAMPLE 4 

14 The mathematical treatment is changed from Multiple Linear Regression to 

15 

Partial Least Squares for calibration transfer to instrument-k from instruments. 

16 

1? In this example, instrument-i is a NIR Systems 6500 bench top instrument 

18 operating in transflectance mode using fiber optics. The instrument is calibrated 

19 as shown in Example 1 using the 72 gasoline samples. The calibration is made 

20 as before modeling (R+M)/2 octane using the second derivative of absorbance 

21 at 1220, 1196 and 1236 nm as independent variables for the multiple regression 
22 

calibration equation. The calibration transfer standards are the same 25 as shown 

23 

24 in Example 5. 

25 Instrument-k is an uncaltbrated NIR Systems OL-6500 operating in 

26 transflectance mode using fiber optics. The 25 calibration transfer standards are 

27 run on this instrument using the calculated octanes, from instrument-!, as the 
dependent variable. This invention makes it possible to calibrate instrument-k 
using a different mathematical approach from that which is used on instrument-i 

31 (multiple linear regression). Partial Least Squares (PLS) is chosen as the 

32 alternative calibration method for instrument-k. The principles' of PLS are well 

33 known (Harold Martens and Tormod Naes, Multivariate Calibration, John Wiley & 

34 Sons; New York, 1989 [ISBN 471-90979-3]). 
35 

PLS is accomplished by modeling the dependent variable of calculated 

36 

37 octane for the 25 calibration standards second derivatives of absorbance over the 

38 wavelength range of 1190-1240 nm (in 2 nm intervals, i.e., 1190, 1192. 1194 . 

39 1240). Cross validation is accomplished using four subset groups. Models are 

40 
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developed using one, two, three, four, five, and six latent variables. The four latent 
variable model seems to provide the best estimate of pump octane for instrument-k 
(i.e., accounted for the most variance without overrating the data), and a 
correlation of 0.9994 with instrumerrt-i's calculated octane vaiues. The far 

latent variable PLS model is used to predict the values of the first 10 samples of 
the original 72 sample calibration set The predicted values show a substantial 
bias of -0.433 (R+M)/2 octane numbers, and a standard error of 0.468. A slope 
adjustment of 0.940 and an intercept adjustment of 5.77 (R+ M)/2 octane numbers 
is made. The remaining 62 samples are predicted as unknowns. The 62 samples 
show a bias of 0.0156 octane numbers and a standard error of prediction of 0.153 
octane numbers for the unknowns. 

it should be noted that a different area of the spectra other than 1 190-1240 
can be chosen, provided that it is functionally related to octane with a reliability at 
least as good as that for measuring such octane directly. This also dearly 
demonstrates that the wavelengths used in calibration equation-1 need not be the 
same in equation 2 which shows some of the utility of the invention. 

EXAMPLE 5 

This example discloses transferring a calibration equation suitable for MTBE 
from one refinery to another without MTBE from one knock engine to another, 
from one instrument after a configuration and lamp change to another instrument, 
with a substantially decreased error of prediction, by means of a double calibration 
transfer. 

This example dramatically demonstrates a utility of the present invention. 
For example, instrument-k is fully calibrated using only gasolines which contained 
no MTBE (methyl t-butyl ether). Recently, Congress passed the reformulated fuels 
act and since instrument-k operates in a refinery that is in a CO (carbon monoxide) 
non-attainment area, the reformulated fuels act specifies that this refinery must 
begin adding MTBE (or some oth r oxygenated feedstock) to its gasoline blends 
at some future date. Because instrument-k is calibrated without the use of MTBE, 
it is expected that substantial errors would result if instrument-k's current 
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calibrstion is used on unknown samples which contain MTBE. To keep from 
completely recalibrating instrument-k, one would want to transfer a calibration 
equation which can accommodate MTBE as it is added to the refineries gasoline 
blending stocks at a future date. In this example, instrument-k is a NIR Systems 
OL-6500 operating in transflectance mode using fiber optics. 

A NIR Systems 6500 is chosen to operate as instrument-! because it was 
once calibrated using 226 samples, some of which contained MTBE, and can 
accommodate samples containing 0-15% by volume MTBE. 

The problem is somewhat compounded by the fact that the only calibration 
standards run on instrument-i during its calibration with MTBE were the eight pure 
compounds (and calculated octanes) shown in Table 6. These eight compounds 
are not run on instrument-k, so a direct calibration transfer is not possible. 

TABLE 6 



I Isooctane 


98.59200 


n-decane 


42.46400 


1 -octane 


57.33100 


2,3,4-trimethylpentane 


100.37100 


n-heptane 


58.10600 


Cyclohexane 


32.51100 


| n-hexane 


66.06500 


I Toluene 


101.69200 



To accomplish the transfer, the eight calibration standards (Table 6) are run 
on instrument-i again. Since the calibration of instrument-i for the MTBE 
calibration, instrument-i has been changed from transmission measurements at a 
path length of 20mm to transflectance measurements at a path length of 18mm 
using fiber optics. The lamp and detector gain are changed as well. 

Under the new set of conditions, using the eight calibration standards, the 
calibration for MTBE is transferred from instrument-i in transmission to instrument-i 
in transflectance. This equation, without slope and bias or other methods of error 
correction (e.g., least squares) is then used to calculate predicted octanes for the 
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25 calibration standards which are also run on instrument-! in transflectance and 

5 instrument-k in transflectance. 

6 Multiple linear regression on the 25 calibration standards Table 7 on 
instrument-k using the calculated octane values as dependent variables and the 
second derivative of absorbance at 1220, 1196, and 1236 nm, yielded 
instrument-k's new calibration equation which could be used to predict the octanes 

1 1 of gasoline blends containing MTBE. 
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The new equation is slope and bias adjusted on ten gasoline samples. Four 
of these samples are premium unleaded gasoline blends [(R+M)/2 octane <* 92] 
and the remaining six are regular unleaded gasoline blends [(R+M)/2 = 87)]. 
These samples show a bias of -0.919 and a standard error of 0.946 (R+M)/2 
octane numbers. 

From this data, a slope correction of 0.950 and a bias correction of 5.30 is 
made. The resulting equation shows a bias of 0.0401 and a standard error of 
prediction of 0.122 on 30 unknown gasoline samples. Some of the pertinent data 
is summarized below. 

SUMMARY TABLE 



Parameter 


Original 
Calibration 
on 

instrument-i 
in 

transmission 


Equation on 
instruments for 
8 

calibration 
transfer 
standards 
in transflectance 


Equations on 
instrument-k for 
25 calibration 
transfer 
standards 
in transflectance 


Final 
calibration 
equation n 
instrument-k 
after slope ami 
bias correction 


A 


87.925 


87.318 


86.238 


87.226 


B, (1220) 


69.864 


94.418 


102.601 


97.471 


B 2 (1196) 


18.581 


23.340 


22.488 


21.364 


B 3 (1236) 


27.644 


28.058 


29.838 


. 28.346 


Correlation 
with 

(R + M)/2 
| octane 


0.992 


* 


0.997 


0.998 


I Standard 
3 Error of 
I Prediction 


0.330 


* 


0.946 


0.122 


| Bias 


0 


* 


-0.919 


0.0401 I 



The Summary Table shows a very interesting phenomenon in that the 
transferred calibration equation is approximately two times more accurate than the 
original equation. This would appear to stem solely from the fact that the octanes 
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3 

4 used in the slope and bias adjustment for instrument-k are more precisely known 

5 than the calibration samples used on instrument-i. 

6 This example also goes against a commonly held opinion of some NIR 
spectroscopists. These scientists believe that an equation developed for one 
refinery is not useful for predictions at another. In this example, the MTBE 
calibration for instrument-i is developed at Ashland's Catlettsburg refinery. 

, , Catlettsburg is a very complex refinery by industry standards and has 12 gasoline 

12 blending components. Instrument-k, on the other hand, analyzes data at 

13 Ashland's St. Paul Park, Minnesota, refinery which has only six gasoline blending 
components and is relatively much less complex. Moreover.the Catlettsburg and 
St. Paul Park refineries have almost no overlap in terms of crude slates. 

17 A schematic flow diagram for this example is shown in Figure^. 

18 ■ 

19 EXAMPLE 6 

This example, schematically outlined in Figure 3, demonstrates that when 
Y M is a linear function of Y u that the calibration transferred equation Is as effective, 
or approximately as effective, as the equation resulting from the direct cafibration 

24 transfer of the previous examples; % 

25 0-©-, where s<pf#dj(itod) = Y M(achi4l) ). Here we will consider the linear function of Y u 

26 where Y k>a = 3^^^ + 25 
In a previous example, an (R+M)/2 octane calibration is shown using 226 

samples and predicts the (R+ M)/2 octanes of 8 calibration standards (isooctane, 
3Q n-decane, l-octene, 2,3,4-trimethylpentane, n-heptane, cychohexane, n-hexane and 

31 toluene). This calibration is developed on an NIR Systems Model 6500 operating 

32 in transmission mode at a path length of 20mm. 

33 Subsequently this instrument was configured for fiber optics using 
transflectance at a path length of 18mm after the lamp and detector gain setting 
are changed. The 8 calibration standards are run on the instrument using three 

37 times the octane value determined from the transmission measurements plus 25 

38 as the value of the dependent variable. The results of multiple linear regression 

39 are: 

40 
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Afc = 286.951 



3 

5 B M220 - 283.252 

6 B M196 = 70.020 

7 B,,^ = 84.185 

8 Values are then predicted for ten gasoline samples. 

g 

^ The values of the ten gasoline samples predicted are corrected by 

n three-equivalent methods. Method one Involves subtracting 25 from the value 

12 predicted for the (R+M)/2 gasoline samples, dividing the result by 3, and 

13 performing a slope and bias adjustment. Method 2 involves a slope and bias 

14 adjustment using the actually measured octane of the sample. Method three 
involves using three times the gasoline samples octane plus 25 for the samples 
used in the slope and bias adjustment and correcting the values predicted for 

18 subsequent unknowns by subtracting 25 and dividing by 3. In each of these three 

19 methods, a bias of 0.0281 and a standard error of prediction of 0.240 (R+M)/2 

20 octane numbers is observed using an unknown set consisting of 62 unknowns. 
It should be noted that the final equation for the direct calibration transfer 

(Via predicted 88 Y k,« actual) shows a somewhat lower error and bias than by either of 

24 the previous methods (bias = 0.0237, SEP = 0.223 (R+M)/2 octane numbers). 

25 y m actual multiplied by a constant with some other constant added to it yielded a 

26 value for the dependent variable for the transfer that had an effect on the error 

27 term (Y M p,^^ - Y ls tcblMS ) which is not removed exactly by either the inverse 

28 mathematical function or slope and bias adjustment 
29 

3q However, if we generate the calibration transfer equation (without slope and 

31 bias adjustment) and multiply the regression constants (A*, B M220 , B k1196 , and 

32 B k 1236 ) by three and add 25 to A k , we can arrive at the exact answer that the direct 

33 method (Y Mpredlcted = Y; tactual ) yields after using a slope and bias adjustment using 

34 the actual (R + M)/2 octanes of the ten samples or by slope and bias adjusting the 
35 

equation using three times the value of octane plus 25 and subtracting 25 from 
37 each predicted octane and dividing by 3. 
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EXAMPLE 7 



1 
2 
3 
4 

5 A calibration equation for an uncalibrated instrument is obtained based upon 

6 the calibration equation established for a reference instrument. Solutions for the 

7 following equations determine: S^'s and A k in the linear calibration equation for 
instrument-k, the uncalibrated instrument, having the form: Y e k = A k +TJ j B j k X jC! k ;and 
which relies on or has recourse to the calibration equation: Y e ' = A'+J^Bj'XjJ; for 

1 1 the calibrated or reference instrument, instrument-i. 

12 Mathematical Criterion I, a classical least squares is used. But, since only 

13 three wavelengths are used, J = 1 to 3; and since 4 = total number of members 
in sample-c, c = 1 to 4. Independent variables, X 1 ^, (where "k" identifies 
instrument-k) have values for spectral features derived from the spectral response 

17 of instrument-k at wavelength-j. Each wavelength-j is within the wavelength range 

18 of from 800 to 2500 nm. Here the spectral features are the second derwafive of 

19 the absorbance determined at three wavelengths-J: 1196, 1220, and 1286 hm. 

20 , , ; 

21 
22 
23 



Mathematical Criterion I 
Ec(Y c '-A k -j;B j k X jc k ) =0; Assume Y.' = Y e k 

24 IcYc^c k -A k I c X 1c k -B 1 k .I^ 1e k X 1 ^B 2 k .2^ k X 1 ^B 3 k •I^ k X 1e k = 0; 

25 kY.V^^-B^EoX,.^^ = 0; and 

27 
28 
29 
30 
31 

32 Even after both changing a bulb which would ordinarily change the spectral 

33 characteristics of an instrument, and also reconfiguring the instrument for fiber 

34 optics operating in transflectance mode at a different path length, it was possible 
to recalibrate such instrument by means of the steps of this invention. ~ 

A calibration is developed on instrument-i for pump octane. The calibration 

38 set consists of 226 samples. These samples are analyzed on instrument-i and a 

39 multiple linear regression is performed modeling (R+M)/2 octane at the second 
40 



Note in Mathematical Criterion I, the form of the calibration equation for the 
reference instrument-i need not be the same as that for the uncalibrated 
instrument-k. 
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derivative of absorbance at 1220, 1196, and 1236 nm. Instrument-i is operating 
in transmission mode using a total path length of 20 mm. The regression 
constants obtained are: 



A' 
B,' 

Ba' 



87.925 
69.864 
18.581 
27.644 

The standard error of estimate for the equation is 0.330. It should be 
pointed out that these results are from two different knock engine operators 
running samples on the same knock engines. These results are considered to be 
good since the ASTM allowable limits for (R+ M)/2 are a standard error of estimate 
of approximately 0.35. These samples were in an approximate box-car distribution 
from 84-94 (R+M)/2 octane. 

Based on the multiple regression calibration equation values are predicted 
for four calibration transfer standards. These standards and their calculated 
(R+M)/2 octanes are shown below. 



Substance 


Calculated Octane 


Isooctane 


98.592 


n-decane 


42.464 


l-octene 


57.331 


2,3,4-trimethylpentane 


100.371 I 



It should also be pointed out that instrument-i is calibrated over a 
four-month period (i.e., samples were run as they came into the lab) and this is not 
an optimal way to do a calibration. Also not optimal was the fact that these 
calibration standards were not run on the instrument promptly after or substantially 
at the same time as the initial calibration, but not until more than one month after 
the calibration was completed. 

Subsequent to this, the lamp on instrument-i burns out and a new lamp is 
inserted. The instrument is also configured for fiber optics operating in 
transflectance mode at a different total path length (~ 18 mm). In addition to the 
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3 

4 configuration mode and path length changes, the detector gain has to be changed 

5 to accommodate the new lamp and lower throughput of the fiber optic 

6 configuration (as opposed to transmission). 
The calibration transfer standards are rerun on instrumeni-i and a 

transferred calibration equation is developed modeling calculated octane as a 
function of the second derivative of absorbance at 1220, 1196, and 1236 nm. 

1 1 The regression constants obtained were: 

12 A k 87.717 

13 B, k = 92.471 

14 B 2 k = 24.079 



/ 
8 
9 
10 
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22 



34 
35 



B„ k - 32.054 



15 

,6 k * 

17 B, is very different from B, 1 whereas A k , B 2 k and B 3 k are within 6 (R+M)/2 

18 octane numbers from A\ B 2 ' and B 3 \ respectively. 

1 9 Subsequently, 97 samples are analyzed for (R + M)/2 octane. None of these 
samples are the same as before. Furthermore they are analyzed by multiple 
operators on different knock engines from those used before. The calibration 

~ transfer equation yields a slight bias of 0.0392 (R+ M)/2 octane numbers (note that 

24 the bias between two different knock engines is commonly much higher than this) 

25 and a standard error of estimate of 0.371 (R + M)/2 octane numbers. A slope and 

26 bias adjustment is made using two samples; one at 87.4 octane and one at 92.5 

27 octane, and yields a slope correction of 0.9402 and intercept correction or bias 

28 

correction of 5.297 and a final equation of: 

29 . 

30 = 87 774 

31 B, k = 86.945 

32 B 2 k = 22.640 

33 B 3 k = 30.139 
The slope and bias adjusted equation has a bias of 0.0 and a standard error 

of prediction of 0.343. 
36 

37 To demonstrate that this could not be accomplished directly from the 

33 original calibration equation from transmission data, the original calibration equation 
39 is used to predict the (R + M)/2 octane for each of the 97 gasoline samples. The 

40 
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result shows a bias of -0.482 octane. A slope and bias adjustment only using ail 



1 

2 -48- 
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4 

5 97 samples on this original calibration equation shows a bias of 0 and a standard 

6 error of estimate of 1.96 (R + M)/2 octane. 

7 Also note that the two point slope and bias adjustment can also be 

8 visualized as a generalized least sum error reducing step in the sense that it can 

g 

be done equivalents using techniques of geometry. For example, defining a line 
1 y extending through the actual values measured, one can determine constants to 

12 multiply the line extending between the values predicted before the adjustment in 

13 such a manner as to make the two lines equivalent. Hence this is generalized 

14 least sum because the approach does not use squares. 

15 
16 

1? MODIFICATIONS 

18 Specific compositions, methods or embodiments discussed are intended to 

19 be only illustrative of the invention disclosed by this specification. Variations on 

20 these compositions, methods or embodiments are readily apparent to a person of 

21 skill in the art based upon the teachings of this specification and are therefore 
22 

intended to be included as part of the invention disclosed herein. 

23 

24 For example, an equivalent method to those disclosed in the examples, 

25 intended to be within the claimed scope of this invention, is determining the bias 

26 value for, for example, A* in calibration equation-k for instrument-k based upon 

27 calibration equation-i for instrument-i, wherein the value for A 1 Is assumed to be 
equal to some arbitrary value, such as zero, or a value determined in connection 
with A* and a function of A 1 ; i.e., a partially calibrated instrument. The "correct" 

31 value for A\ is that value which when used in a calibration equation-k determines 

32 the value for a dependent variable-k ( which is substantially equal to the actual value 

33 of a physical or chemical property desired. The "correct" value for A k can be 

34 obtained from employing a slope and bias correction to that equation which 
resulted from employing an arbitrary value for A' or a value determined by some 
function of A*. Another equivalent method like above occurs when the final 

38 equation that would otherwise result from a slope and bias equation may be 

39 disguised by using a function of either m or b, (the terms discussed in Example 6). 



28 
29 
30 



WO 94/08225 PCT/US92/08547 



13 
14 
15 
16 



20 
21 
22 



27 
28 
29 



34 
35 
36 



-49- 



1 
2 
3 

4 This will clearly change the calculated or predicted value for the dependent 

5 variable. However, just as was shown in the specification that a function of Ys 

6 from calibration equation-i can be used, so also can a function of m or 6 be used 

7 and then the resulting dependent variable transformed back to a vaiue that wiii be 
substantially equal to the actual value for a chemical or physical property of 
interest. 

1 1 Still another variation within the scope of this invention involves how values 

12 are assigned to dependent variable-k. These values for dependent variable-k, that 
can be used to determine constants-k in calibration equation-k or for performing 
a slope and bias correction, can be determined without using a primary test 
method but instead by reference to calibration equation-i and spectral data-i 

17 obtained from instrument-i .or possibly any other calibrated instrument in 
i s conjunction with its corresponding calibration equation which involves substantially 
19 the same range of interest and relevant spectral data. 

Still another variation of this invention is to use different spectral features in 
calibration equation-i from those in calibration equation-k, e.g., different 

23 wavelengths or different mathematical transforms or both; provided such 

24 differences are each correlatable to a physical or chemical property, of interest. 

25 Reference to documents made in the specification is intended to result in 

26 such patents or literature being expressly incorporated herein by reference, 
including any patents or other literature references cited within such documents. 

In most examples, only three wavelengths were used. The number, 
3Q however, can vary depending upon the statistical reliability sought or required. If 

31 the indirect method has an SEP of v for some property which is smaller than the 

32 interval of reliability for the direct method for determining such property, then there 

33 is a practical limit to how small it is reasonable to seek to make V, and how 
many wavelengths should be required. 

Another modification relates to the spectral features which can be derived 
3? from absorption spectra. In the examples, both second and higher derivatives of 

38 intensities versus wavelengths or absorbances, were used. Clearly, higher 

39 derivatives such as first, second, third, fourth, and fifth or higher are also expressly 

40 
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4 contemplated for this invention. Values for a spectral feature can equal simply the 

5 value found for a particular wavelength, or a range of wavelengths around some 

6 reference wavelength; e.g., absorbances. Spectral features are not intended to 

7 suggest that the spectral data is necessarily different in value from a spectral 
g 

feature. It is clearly contemplated that it may be more complicated in that a linear 

9 

combination of spectral features or mathematical transforms thereof are within the 

10 

1 1 intended scope of this invention. 

12 Also contemplated as a variation of this invention is a linear combination of 

13 still higher adjacent order derivatives, such as third and fourth order derivatives or 

14 higher adjacent order derivatives, of spectral data. See Examples for an 
example. 

16 

17 Further, any set of members or range of numbers cited in the specification 

^8 is intended to incorporate expressly herein any member or number falling within 

19 such set or range, respectively, including any subset or range within any set or 

20 range cited. 

21 What is claimed is: 
22 

23 

24 

25 

26 

27 

28 
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31 
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CLAIMS 



1 
2 
3 
4 

5 1. A method for calibrating a second instrument by determining 

6 constants in a calibration equation, calibration equation-2, for said second 

7 instrument by recourse to a first instrument having a calibration equation, 
8 

calibration equation-1 , comprising in any order: 

9 
10 

1 1 (a) Determining from calibration equation-1 a value indicative of each 

12 dependent variable- 1 corresponding to each member of a calibration set based on 

13 spectral data-1 for each member of said calibration set measured by means of 

14 instrument- 1 ; and 

15 
16 

1 7 (b) Determining constants-2 in said calibration equation-2 so that the 

18 sum of absolute differences between: 
19 

20 (1) said value indicative of each said dependent variable-2 for each member 

2 1 of said calibration set based upon calibration equation-2 and data-2 corresponding 



to each said member; and 



is substantially minimized. 



22 
23 
24 

25 (2) that corresponding value of dependent variable- 1 calculated for each 

26 said member based upon calibration equation-1 and data-1 corresponding to each 

27 said member; 

28 
29 
30 
31 

32 2. The method of claim 1, wherein said calibration set has at least one 

33 member that yields spectral data-2 when measured by instrument-2 which in 

conjunction with calibration equation-2 determines a value for said dependent 

variable-2 for said at least one member which is indicative of a value which is 
36 

^ significantly different from that value determined by a direct test method for which 
33 said calibration equation-1 and calibration equation-2 were established to model. 

39 
40 
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3. The method of Claim 1, wherein calibration equation- 1 has the 

5 following form: 

6 Y 1 8 = A 1 + Ej (B'j • I'jJ 

7 wherein: 
8 



9 
10 



Y 1 , is dependent variable-1 for sample member s; 
A 1 is an intercept constant- 1 

11 B 1 j is a slope constant-1 for every sample member s at wavelength 

12 j; 

13 l 1 j8 is an independent variable determined from spectral data-1 

14 corresponding to sample member s at wavelength j; 
15 

16 

17 and said minimization of said absolute differences is achieved by substantially 
<ig minimizing the following mathematical expression: 

19 

20 Ls^-A'-S^.l^} 2 . 
21 

22 4. The method of Claim 3, wherein said independent variable is 

23 

determined from spectral data-1 by taking a mathematical transform of 
absorbance. 

25 
26 

27 5. The method of Claim 4, wherein said mathematical transform 

28 comprises a differential of an absorbance spectra 
29 

30 
31 
32 
33 

34 Y 1 , = A 1 + £ (B 1 , • g 

35 wherein: 

Y 1 S is dependent variable-1 for sample member s; 
A 1 is an intercept constant-1 

39 B 1 - is a slope constant-1 for every sample member s at wavelength 

40 i' 



36 
37 
38 



6. The method of Claim 1, wherein said calibration equation-1 and 
calibration equation-2 have the same form shown hereinafter: 
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Y Z . - A 2 + 5 (B 2 j • ly 



wherein: 



i 

2 
3 

4 l' js is an independent variable determined from spectral data-1 

5 corresponding to sample member s at wavelength j; and 
6 

7 
8 

9 " 

10 Y 8 is dependent variable-2 for sample member s; 

11 A 2 is an intercept constant-2 

12 B 2 is a slope constant-2 for every sample member s at wavelength 

13 j; 

14 i2 • 

^ I j S is an independent variable determined from spectral data-2 

16 corres Po n d'ng to sample member s at wavelength j; and said determining said 

17 minimization of said sum of absolute differences comprises setting slope and 
1 s intercept constants for calibration equation- 1 to values established by the following: 

19 

20 B'j = B 2 { £, (l 2 j9 • IV > / { L (" 1 js • 'V> > and 

21 
22 
23 
24 
25 

26 7 - Tn® method of Claim 1, wherein said calibration equation- 1 and 

27 calibration equation-2 have the same form shown hereinafter: 
28 
29 
30 
31 
32 

33 Y 1 6 is dependent variable-1 for sample member s; 

34 A 1 is an intercept constant-1 

35 B 1 j is a slope constant-1 for every sample member s at wavelength 

36 

j> 

37 , 

I j8 is an independent variable determined from spectraJ data-1 
3g corresponding to sample member s at wavelength j; and 

40 



A 1 = X. { Y 2 , - 5 (B 1 . • } / I 

where I is the number of members in said calibration set. 



Y 1 S = A 1 + £j (B 1 , • IV. 

wherein: 
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4 Y% = A 2 + £ (B* • \\) 

5 wherein: 

6 Y 2 6 is dependent variable-2 for sample member s; 

7 A 2 is an intercept constant-2 

B 2 is a slope constant-2 for every sample member s at wavelength 

d 

.10 j: 

n l 2 js is an independent variable determined from spectral data-2 

12 corresponding to sample member s at wavelength j; and said determining said 

13 minimization of said sum of absolute differences comprises setting slope and 

1 4 intercept constants for calibration equation-1 to values established by the following: 

]* 7 B\ = B 2 { £. I 2 js / £ IV, } and 

18 A 1 = £ { Y 2 S - £ B 1 , • I 1 ,, > /n 

19 

20 where n is the number of members in said calibration set. 
21 

22 

23 8. The method of Claim X, wherein there is a further improvement of 

24 applying a slope and bias correction comprising the step of substituting values for 

25 m and b into the following expression: 

26 
27 
28 

29 ' 

3q wherein m and b are defined by the following expressions: 

31 

32 m = {(I.Y, 8 JXVn-LO^ Y e k )}/{{^Y c k E^/n-E^" Y e k } 

33 



Y\ = { m x A(1) + b } + £ { m x B\ x } 



34 b = i/ndcY^-m^Y,"} 



where: 

36 



37 Dependent-variable-(ac), Y%, is selected from the group of values consisting 

38 of: a directly determined value for property-p of a member c; a predicted value for 

39 property-p of member c, by means of a calibration equation for property-p within 
40 



» 
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a range said calibration equation is established to model; wherein each value for 
said property-p is in a range that said calibration equation-1 was established to 
model; and a combination thereof; 
and 

n is the number of members to a set of slope and bias samples that yield 
calculated values indicative of values within a range that said calibration equation-1 
is intended to model. 

10. The method of Claim 1, wherein said spectral data-1 comprises 
absorbances or mathematical transforms thereof. 

11. The method of Claim 10, wherein said dependent variable-1 
corresponds to a value indicative of octane or cetane and said absorbances or 
mathematical transforms thereof involve frequencies in the range of 800 to 2500 
nanometers. 

12. The method of Claim 11, wherein at least one member of said 
calibration set is selected from compounds active in the range of 800 to 2500 
nanometers. 

13. The method of Claim 12, wherein said compounds active in the range 
of 800 to 2500 nanometers includes materials selected from the group consisting 
of: n-paraffins, iso-paraffins, aromatics, alcohols, olefins, halogenated alkanes, 
ketones, esters, ethers, and naphthenes. 

14. A method for calibrating or recalibrating a first instrument, instrument- 
k, by reference to a second instrument, instrument-i, comprising in any order: 

(a) obtaining spectral data by means of instrument-k and instrument-i, 
corresponding respectively to data-k and data-i for each member of 
a calibration set; 



WO 94/08225 



PCI7US92/08547 



-56- 

(b) determining by means of a calibration equation, calibration equation- 



(c) determine a value, dependent variabie-f, related by function-f to said 



i 

2 
3 
4 

5 i, of said instrument-i, a value, value-i, for a dependent variable, 

6 dependent variable-i, corresponding to each member of said 

7 calibration set; 

8 
9 
10 

n value-i for said dependent variable-i, for each member of the 

12 calibration set; 

13 

14 (d) selecting values for said constants-k in calibration equation-k so that 
15 

the sum of absolute differences between: 

16 
17 

13 (1) each dependent variable-k calculated based on said 

19 calibration-k with said selected constants-k and data-k for 

20 each member of said calibration set; and 
21 

22 

23 

24 said calibration set determined by function-f and 

25 corresponding value-i for dependent variable-i of each said 

26 member; 

27 is substantially minimized. 

28 . 

29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 



(2) corresponding dependent variable-f for identical members of 
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(b) Determining constants-2 in said calibration equation-2 so as to 
minimize the sum of absolute differences between: 
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2 . AMENDED CLAIMS 

[received by the International Bureau on 10 August 1993 (10.08.93); 

3 original claims 1-14 replaced by amended claims 1-10 (5 pages)] 
4 

5 1 . A method for calibrating a spectral instrument to determine a value 

6 of a property of an unknown sample by using constants in a calibration equation-2, 

7 for said instrument, by reference to a spectral instrument having a calibration 

8 equation- 1, said method comprising in any order: 
9 

10 

1 1 (a) Determining from calibration equation- 1 a value of a dependent 

1 2 variable- 1 corresponding to each member of a set of calibration transfer standards 

13 based on spectral data-1 for each member of said set of calibration transfer 

14 standards; and 
15 

16 
17 
18 
19 

20 (1) Each said value of dependent variabie-2 for each member 

21 of said set of calibration transfer standards using spectral data-2 

22 

measured by said instrument; and 

23 

24 '" 

25 (2) That corresponding value of dependent variable-1 

26 determined for each member of said set of calibration- transfer 

27 standards determined by calibration equation- 1. 

28 
29 
30 
31 
32 

33 2. A method for calibrating a second spectral instrument to determine 

34 a value of a property of an unknown sample by using constants in a calibration 

35 equation-2, for said second instrument, by reference to a first spectral instrument 

36 

having a calibration equation- 1, said method comprising in any order: 

37 

38 
39 
40 



(c) Determining values of said property in said unknown sample by 
use of said instrument 
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(2) That corresponding value of dependent variable-1 



(a) Determining from calibration equation- 1 a value of a dependent 
variable-1 corresponding to each member of a set of calibration transfer standards 

5 

6 based on spectral data-1 for each member of said set of calibration transfer 

7 standards measured by means of instrument-1 ; and . 
8 

9 (b) Determining constants-2 In said calibration equation-2 so as to 

10 minimize the sum of absolute differences between: 
11 

12 

13 (1) Each said value of dependent variable-2 for each member 

14 of said set of calibration transfer standards using spectral data-2 

15 measured by said second instrument; and 
16 

17 

18 

19 determined for each member of said set of calibration transfer 

20 standards measured by said first instrument. 
21 

22 (c) Determining values of said property in said unknown sample by 

23 

use of said second instrument 

24 
25 

26 3. A method for calibrating or recalibrating a first instrument, instrument- 

27 1, by reference to a second instrument, instrument-2, over a range of valued for 

28 a property of Interest comprising in any order: 

29 
30 

(a) obtaining spectral data by means of instrument-1 and instrument-2, 

32 corresponding respectively to data-1 and data-2 for each member of 

33 a set of calibration transfer standards; 

34 

35 (b) determining by means of a calibration equation, calibration equation- 

36 1 , of said instrument-1 , a value for a dependent variable, dependent 

37 

variable-1, corresponding to each member of said set of calibration 
39 transfer standards; 
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3 

4 (c) determining a value for dependent variable-2 that is relatable by a 

5 function-f to said property of interest, for each member of the set of 

6 calibration transfer standards; 



28 
29 
30 
31 



35 
36 
37 



(d) selecting values for said constants-1 in calibration equation-1 so that 
the sum of absolute differences between: 



8 

9 

10 

ii 

12 0 ) each dependent varlable-1 estimated by said calibration-1 with 

13 said selected constants-1 and data-1 for each member of said 

14 set of calibration transfer standards; and 

15 
16 

17 W corresponding dependent variable-2 for identical members of 

18 s©* ° f calibration transfer standards, 

19 

20 is substantially minimized. 
21 

^ 4 - A method according to Claim 1, 2 or 3, wherein said set of calibration 

24 transfer standards have values of said dependent variable which are different from 

25 the value of said dependent variable which would be obtained by a primary 

26 reference test method, and wherein said set of calibration transfer standards have 

27 compositions different from the composition of said sample, and wherein the 
determining of constants-2 is by a method which comprises: 

(a) partial least squares; 

(b) multiple linear regression; and/or 
32 (°) principal component regression. 
33 

34 5. A method according to Claim 1 , 2 or 3, wherein said spectral data-1 and 

the spectral data-2 comprise infrared spectral data, and wherein said sample 
comprises liquid, and wherein said sample comprises liquid, and wherein said set 

3 8 calibration transfer standards comprises at least one composition different from 

39 those used to calibrate instrument- 1. 

40 
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8. A method according to Claim 1, wherein calibration equation-2 has 



6. A method according to Claim 1, wherein said set of calibration transfer 
standards comprises at least one standard which is a pure compound or mixture 
6 of pure compounds. 

7 

8 7. A method according to Claim 1, wherein said set of calibration 

9 transfer standards have at least one member that yields spectral data-2 when 
™ measured by instrument-2 which in conjunction with calibration equation-2 

12 determines a value for said dependent variables for said at least one member 

13 which is incficative of a value which is significantly different from that value 

14 determined by a direct test method for which said calibration equation- 1 and 

1 5 calibration equation-2 were established to model. 
16 

17 

18 

ig the following form: 

20 Y 2 t = A 2 + 5 (B*j ly 

21 wherein: 

22 Y 1 s is dependent variable- 1 for sample member s, based upon 
23 

calibration equation- 1 for instrument-1; 

24 

Y z 8 is dependent variable-2 for sample member s; 

25 

26 A 2 is an intercept constant-2 

27 B 2 is a slope constant-2 for every sample member s at wavelength 

28 i; 

29 l 2 j8 is an independent variable determined from spectral data-2 
30 

corresponding to sample member s at wavelength j; 

32 

33 and said minimization of said absolute differences is achieved by substantially 

34 minimizing the following mathematical expression: 

35 

36 {Y\ - A 2 - £j (B 2 • l 2 ^ } 2 ; and wherein said independent variable 

37 is determined from spectral data-2 by taking a mathematical transform of 

38 absorbance, and wherein said mathematical transform comprises a differential of 
an absorbance spectra. 

40 



22 
23 
24 



m = {(^VlY c Vn-L(Yc*Y c k )>/{{I c Y c k I c Y c h }/n-^Y e fc Y c k } 
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4 

5 9. A method according to any preceding claim, wherein there is a further 

6 improvement of applying a slope and bias correction, whereby the accuracy of 

7 dependent variable, Y 1 ., is increased to Y 1 .. comorisina the sten nf 

8 values for m and b into the following expression: 
9 

J° y\ - { m x A(1) + b } + £ { m x B 1 , x 1^ > 

12 

13 wherein m and b are defined by the following expressions: 

14 
15 
16 

]l b = 1/n{E 0 Y c '-mI 0 Y e ,l > 

19 where: 

20 Dependent-variable-(ac), Y* c , is selected from the group of values consisting 

21 of: a directly determined value for property -p of a member c; a predicted value for 
property-p of member c, by means of a calibration equation for property-p within 
a range said calibration equation is established to model; wherein each value for 

25 said property-p is in a range that said calibration equation-1 was established to 

26 model; and a combination thereof; 

27 and 

28 n is the number of members to a set of slope and bias samples that yield 

29 

calculated values indicative of values within a range that said calibration equation-1 

30 

mm is intended to model. 
31 

32 

33 10. A method according to Claim 1 through 8, wherein there is a further 

34 improvement comprising an error reducing step. 
35 
36 
37 
38 
39 
40 



WO 94/08225 



PCT/US92/08547 



1/3 



eox i 

CALBRATE NSTRUMENT-1 OR i 
BY OEFttC A CMJBRATKW 
EQUATON-I WHJCH MODELS 
THE CORRESPONDENCE BETWEEN 
SPECTRAL DATA MEASURED 
BY INSTRUMENT-! AM) PI 
FOR EACH MEMBER OF THE 
SET OF SAMPLES WITH AN SEE 
NOT SGNFCANTLY 
GREATER THAN THAT SEE 
►KRENT N THE DFECT TEST 1 
METHOO USED TO OETERUNE P 1. 


FIG. 1 


BOX 2 

DEJTRMJC SOME PROPERTY 
Of WTEREST ( P 1) *ITH A 
DIRECT TEST METHOO FOR 
EACH MEMBER OF THE SET OF 
SAMPLES WITH A STANDARD 
ERROR OF ESTMATE. SEE. 
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BOX 3 

DETERUK A DEPENDENT 
VARIABLE VALUE 
(DEPENDENT VARUSLE-1) 
FOR EACH MEMBER OF A 
CALIBRATION SET BASED 
UPON ITS SPECTRA MEASURED 
BY MSTRUMENT-1 AND 
CALBRADON EOUATION-1 



T 



BOX 5 

SELECT CONSTANTS ( CONSTANTS-2 
OR k IFOR CALBRADON EOUATIONS-2 
SO THAT THE SUM OF ABSOLUTE 
OFFERENCES FOR EACH MEMBER 
OF SAP CALBRATPM SET BETWEEN 
CORRESPONDMG DEPENDENT 
VARWBLE't AND DEPENDENT 
VARIABLE-2 ARE MMMZED: E.G. 
SUM OF ALL I DEPENDENT 
VMABLE-1 MMUS DEPENDENT 
VARBLE-2 iZ FOR EACH MEMBER 
OF THE CALSRATION SET 6 
MNUZED. 



BOX 4 
OBTMN SPECTRAL OATA FOR 
EACH MEMBER OF SAD 
CALBRATUN SET BY MEANS 
OF MSTRUCNT-2 AND 
DK3TIZE OR MAKE AN ANALOG 
SBNAL OF SUCH OATA TO OETERUNE 
A SET OF MDEPENDENT 
VARIABLES ( MEPENDENT 
VAWABLES-2 > 



BOX 7 

PREDICT P I. FOR • UNKNOWNS 
WITH A STANDARD ERROR OF 
PREDOOIKSEP) THAT IS NOT 
SKSNfflCANTLY GREATER THAN 
•SEE " 



BOX 6 

ERROR REDUCE, to. 
SLOPE AND BUS ADJUST. 
CONSTANTS-2 WITH KNOWNS 
ON NSTRUMENT-2 



BOX 8 
REPEAT FOR SAME OR 
OFFERENT INSTRUMENT 
F OEStfSD 



BOX 9 
FINISHED 
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EXAMPLES 
FLOW DIAGRAM 



Instrument-i 

calibrated in transmission mods 
capable of prodding MTBE blends 



® 



lnstrument-4 in transflectance mode 
calibrated by transfer of calibration 
equation from® 



I 



© 



Instrument-* in transflectance mode 
calibrated by transfer from® 
capable of predicting MTBE blends 



I 



lnstrwmnt-4 In Cttkmbwg RHIimy 
Uwtnumnt-k In St Paul Ratlnary 



* 



Figure 2 
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EXAMPLE 6 
Flow Diagram 
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